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Abstract 

This  research  models  and  analyzes  methods  to  damage  Bacillus  spores  through 
various  heat-treatments.  AFIT  researchers  have  been  examining  and  developing  meth¬ 
ods  to  characterize  the  effects  of  heating  Bacillus  anthracis  spores  to  high  tempera¬ 
tures  and  for  short  durations.  The  current  laboratory  experiment  being  conducted  is 
designed  to  imitate  the  neutralization  of  Bacillus  spores  by  a  thermal  pulse  similar  to 
that  of  a  conventional  weapon  detonation.  The  research  contained  in  this  report  de¬ 
signs  a  thermal  model  that  replicates  the  current  laboratory  experiment  and  evaluates 
the  rate  of  thermal  diffusion  throughout  the  Bacillus  spores. 

In  addition,  a  micro-model  of  the  effects  of  dry  and  wet  heating  on  a  spore 
is  presented.  By  applying  heat  to  a  spore,  we  energize  adsorbed,  absorbed,  and 
chemically  bound  water  molecules.  These  energized  molecules  have  greater  mobility 
and  diffusion  within  the  spore,  as  well  as  between  the  spore  and  the  surrounding 
environment.  The  water  release  permits  hydrolysis  reactions  to  take  place  with  the 
spore’s  DNA  and  other  proteins.  These  chemical  reactions  degrade  the  DNA  and 
proteins  to  such  an  extent  that  the  DNA  cannot  be  repaired  or  replicated,  thus  causing 
spore  death.  We  further  assert  that  damage  to  a  spore  is  based  on  a  certain  initial 
DNA  information  content  and  the  spore  population’s  protein  ‘fitness’.  Once  this 
protein  fitness  level  is  degraded  below  a  critical  threshold,  the  protein  can  no  longer 
repair  DNA  and  thus  the  spore  is  unable  to  germinate,  or  produce  outgrowth.  These 
assertions  allow  us  to  create  a  probability  of  kill  model  based  on  water  mobility, 
hydrolysis,  a  spore’s  DNA  information  content,  and  the  spore  population’s  protein 
fitness. 
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MODELING  THERMAL  INACTIVATION  OF  BACILLUS 

SPORES 


I.  Introduction 


1.1  Background, 

In  1979,  in  Sverdlovsk,  Russia,  a  Soviet  bioweapons  facility  unintentionally  re¬ 
leased  spores  of  Bacillus  anthracis  ( B.a .),  the  bacterium  that  causes  anthrax.  This 
led  to  an  epidemic  of  inhalational  anthrax  in  Russia  and  18  cases  of  exposure  to  an¬ 
thrax  in  the  United  States  (U.S.).  In  1995,  a  Japanese  cult,  Aum  Shinrikyo,  released 
aerosols  of  anthrax  and  botulism  throughout  Tokyo  and  within  the  Tokyo  subway 
system  at  least  eight  times.  Although  the  strain  of  dispersed  anthrax  was  used  for 
animal  vaccinations  and  was  not  significantly  dangerous  to  humans,  the  attacks  pro¬ 
duced  widespread  panic.  In  2001,  following  the  terrorist  attacks  of  11  September, 
B.a.  spores  were  mailed  through  the  U.S.  Postal  Service  to  members  of  the  U.S. 
Senate  and  to  national  news  agencies.  These  biological  attacks  resulted  in  22  com¬ 
bined  inhalational  and  cutaneous  cases  of  anthrax  infection.  Of  these  22  cases,  five 
Americans  died  and  “the  nation  was  terrorized  in  what  became  the  worst  biological 
attacks  in  U.S  history”  [3].  These  attacks  illustrated  the  prevalent  need  of  the  U.S. 
to  consider  the  detection,  recognition,  environmental  surveillance,  decontamination, 
and  destruction  of  anthrax. 

The  Working  Group  on  Civilian  Biodefense  identified  anthrax  as  a  biological 
agent  that  may  be  used  as  a  weapon  and  “in  worst  case  scenarios,  could  cause  disease 
and  deaths  in  sufficient  numbers  to  gravely  impact  a  city  or  region”.  At  least  13 
countries  are  suspected  to  contain  offensive  biological  weapons  programs  [4],  There¬ 
fore,  both  political  and  military  factions  within  the  U.S.  have  stressed  the  need  for 
counter  terrorism  weapons  and  biodefensive  systems  capable  of  destroying  biological 
and  chemical  weapons  [47]. 
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1.2  Purpose  of  the  Research 

Although  multiple  data  sources  can  be  found  regarding  dry  heat  damage  to 
spores  over  long  periods  of  time,  very  little  research  has  been  completed  on  extremely 
short  duration  exposures  consistent  with  detonation  of  a  conventional  weapon.  Also 
extremely  important  is  understanding  the  effectiveness  of  weapons  against  biological 
targets.  Researchers  at  the  Air  Force  Institute  of  Technology  (AFIT),  working  in 
conjunction  with  the  Air  Force  Nuclear  Weapons  and  Counterproliferation  Agency 
(AFNWCA),  has  been  researching  and  developing  methods  to  characterize  the  effects 
of  heating  Bacillus  anthracis  spores  to  high  temperatures  and  for  short  durations. 
The  current  laboratory  experiment  being  conducted  by  AFIT  researchers  is  designed 
to  replicate  the  neutralization  of  B.a.  spores  by  a  thermal  pulse  similar  to  that  of 
weapon  detonation.  Specifically,  the  objective  of  the  experiment  is  to  determine  the 
probability  of  kill  for  B.a.  spores  exposed  to  high  temperatures,  300  to  1300  °K  (27  to 
1027  °C),  for  short  duration  exposures,  0.01  to  3  seconds  [5].  The  experiment  consists 
of  rastering  a  solid  state  laser  (NckYAG  laser)  beam  to  indirectly  heat  the  spores  for 
various  time  intervals  and  temperature  ranges  using  a  laboratory  produced,  custom 
black  body  radiator.  The  spores  are  analyzed  before  and  after  the  laser  heating  to 
determine  the  probability  of  kill  dependent  on  regeneration  of  each  spore.  In  order 
to  accurately  predict  the  effectiveness  of  the  laser  heating  method,  it  is  necessary  to 
understand  how  the  thermal  pulse  (generated  by  the  laser)  propagates  through  the 
test  apparatus.  To  accomplish  this,  a  thermal  model  will  be  designed  to  replicate  the 
current  laboratory  experiment  and  evaluate  the  rate  of  thermal  diffusion  throughout 
the  B.a.  spores.  This  model  will  be  numerically  evaluated. 

Once  the  thermal  exposure  is  understood,  it  will  be  necessary  to  characterize 
the  influence  this  has  on  the  spore  itself.  There  have  been  many  proposed  models 
regarding  heat  damage  to  spores  which  are  briefly  summarized  here.  Setlow  theo¬ 
rizes  that  DNA  damage  is  the  primary  cause  of  thermal  damage  because  the  target 
for  dry  heat  treatment  “is  at  least  in  part  DNA,  as  there  is. . .  much  DNA  damage 
and  loss  of  DNA  repair  capacity”  [39].  Jagannath,  Tsuchido,  and  Membre  speculate 
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that  heating  spores  that  are  suspended  in  certain  foods  is  more  efficient  due  to  the 
existence  of  acidic  pH  [24],  Geeraerd,  Herremans,  and  Van  fmpe  hypothesize  that 
the  components  within  the  spores  experience  a  log  linear  inactivation  and  the  spores 
themselves  encounter  a  log10  rate  of  thermal  decay.  Casolari  assumes  a  spore  is  killed 
by  a  “lethal  hit  of  a  water  molecule”  due  to  the  high  amount  of  energy  that  water 
molecules  contain  [18].  One  contention  that  all  these  models  have  in  common  is  that 
spores  are  extremely  resistant  to  both  wet  and  dry  heat.  They  also  agree  that  while 
dry  heat  damage  is  caused  by  at  least  some  DNA  damage,  wet  heat  does  not  produce 
any  DNA  damage. 

We  have  discussed  models  which  theorize  damage  to  Bacillus  spores  by  various 
methods.  These  models  use  multiple  Bacillus  species  such  as  anthracis,  cereus,  and 
subtilis.  Therefore,  it  is  necessary  to  discuss  the  validity  of  using  different  species  to 
characterize  spore  damage  due  to  heat  treatment.  One  of  the  obvious  concerns  in 
regards  to  heat  damage  is  varying  spore  heat  resistance  between  species  but  “these 
mechanisms  appear  to  be  conserved  across  Bacillus  species”  [40].  All  species  of  the 
Bacillus  genus  are  genetically  similar  with  the  exception  of  the  existence  of  an  outer 
spore  layer  called  an  exosporium.  But  the  exosporium  is  not  believed  to  provide 
increased  protection  for  the  spore  [21].  For  these  reasons,  we  assume  that  spore  kill 
mechanisms  are  similar  across  the  Bacillus  genus. 

The  approach  of  this  research  focuses  on  hydrolysis  reactions  within  the  spore. 
Hydrolysis  is  a  reaction  of  a  molecule  with  water.  In  this  case,  water  can  react  with 
biopolymers  such  as  DNA  and  proteins  resulting  in  depolymerization  (breaking  a 
large  strand  of  polymers  into  two  smaller  ones)  or  removal  of  a  side  group  of  the 
polymer  [27].  We  propose  that  these  hydrolysis  reactions  decimate  the  levels  of  DNA 
and  proteins  to  such  an  extent  that  the  DNA  can  no  longer  be  repaired.  This,  in 
turn,  prevents  the  spore  from  germination  and  effectively  kills  the  spore.  In  order 
to  assess  the  probability  of  kill,  a  water  production- diffusion  micro  model  based  on 
hydrolysis  reactions  is  needed  to  determine  thermal  damage  to  Bacillus  spores.  The 
hydrolysis  reactions  occur  mainly  in  the  cortex  and  the  core  of  each  spore  and  since 
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the  water  molecules  are  assumed  to  be  mobile,  this  model  will  focus  on  diffusion  of 
water  throughout  each  section  of  the  spore  in  relation  to  ambient  water  concentration. 
The  ambient  water  concentration  is  dependent  on  whether  the  spore  is  being  heated 
via  wet  or  dry  heat.  A  wet  heating  process  is  one  in  which  the  spore  is  suspended  in 
a  saturated  aqueous  environment  while  a  dry  heating  process  applies  hot  air,  which 
contains  little  to  no  water  vapor,  to  the  spore.  Either  heating  process  will  energize 
water  molecules  contained  in  a  spore  and  these  energized  molecules  become  mobile. 
The  consequence  of  an  increase  in  core  water  content  is  “more  rapid  heat  inactivation 
of  spore  core  proteins”  [11],  Our  model  will  consider  the  production-diffusion  of  water 
throughout  the  spore  and  relate  it  to  the  rate  of  hydrolysis  reactions  with  DNA  and 
proteins.  Finally,  a  probability  of  kill  model  that  relates  the  damage  of  the  initial 
DNA  information  content  to  the  decimation  of  the  spore  population’s  protein  ‘fitness’ 
will  be  presented. 

1.3  Research  Objectives 

•  Numerically  model  the  current  AF1T  experiment  of  heating  B.a.  spores  and  the 
rate  of  thermal  diffusion  throughout  the  spores  in  order  to  calculate  the  actual 
temperature  and  durations  experienced 

•  Analytically  model  and  analyze  the  micro-level  damage  within  a  B.  a.  spore  due 
to  production,  diffusion,  and  hydrolysis  reactions  of  water  molecules  throughout 
a  spore. 

•  Relate  the  thermal  model  to  production/diffusion/hydrolysis  reaction  model  to 
develop  a  probability  of  kill  model  based  on  the  threshold  killing  mechanisms 
of  B.a.  spores 

1.4  Overview 

The  objectives  previously  defined  will  be  presented  in  the  following  manner.  The 
development,  numerical  methods,  and  results  of  the  thermal  model  will  be  described 


1-4 


in  Chapter  II.  Chapter  III  explains  the  release  of  water  molecules  within  the  spore  due 
to  heat  application.  It  illustrates  the  diffusion  of  these  water  molecules  throughout 
the  spore  which  promotes  hydrolysis  reactions  to  occur.  Chapter  IV  discusses  our 
proposal  that  spore  damage  is  caused  by  hydrolysis  reactions  since  hydrolysis  will 
damage  a  spore’s  DNA  information  content  and  degrade  repair  proteins  to  such  an 
extent  that  the  DNA  cannot  be  repaired.  It  will  describe  and  present  results  from  the 
kill  probability  model.  The  conclusion,  as  well  as  any  suggested  future  work,  will  be 
contained  in  Chapter  V. 


1-5 


II.  Thermal  Model  and  Analysis 


2.1  Current  Experiment 

In  order  to  analyze  the  consequences  of  high  temperatures  for  short  durations 
on  Bacillus  anthracis(B.a.),  researchers  at  AFIT  have  developed  multiple  laboratory 
experiments  consistent  with  detonation  of  a  conventional  weapon.  This  research  will 
numerically  model  the  heating  and  thermal  diffusion  created  by  the  current  exper¬ 
iment.  To  understand  the  thermal  model  developed  by  this  research,  a  thorough 
description  of  the  experiment  is  required. 

The  current  experiment  was  created  to  heat  B.  a.  spores  using  a  solid  state, 
Nd:YAG  laser.  A  soda  lime  glass  coverslip  was  micro-etched  to  generate  wells  in  the 
coverslip.  The  circular  wells  were  50  microns  (/rm)  wide,  between  3  and  5  /irn  deep, 
and  each  coverslip  contained  a  grid  of  532  total  wells  (38  columns  by  14  rows).  See 
Figure  2.1  for  a  top-level  view  of  the  etched  coverslip  [20]. 


The  wells  were  then  populated  with  the  B.  a.  spores.  The  diameter  of  a  B.  a.  spore 
is  approximately  1  /irn  [38].  Viable  wells  contained  between  1  and  15  spores  each. 
A  soda  lime  glass  slide  was  also  micro-etched  with  wells.  Carbon  black  was  placed 
in  these  wells  as  a  optical  absorber.  Because  the  spores  could  not  be  directly  heated 
via  the  laser,  the  carbon  black  served  as  a  ‘hot  plate’  converting  optical  energy  to 
thermal  energy  and  conducting  heat  into  the  spore  wells.  The  carbon  black  wells 
were  approximately  75  /irn  wide  to  ensure  complete  coverage  of  the  wells  containing 


2-1 


the  spores  (50  /irn  wide).  There  are  1500  /mi  between  each  carbon  black  well  in  each 
column,  row,  and  from  all  edges  of  the  coverslip.  The  glass  slide  was  then  placed  over 
the  etched  coverslip  holding  the  spores.  The  slide  with  the  carbon  black  wells  was 
removed  after  heat  treatment  so  the  spores  could  be  viewed  under  the  microscope. 
See  Figure  2.2  for  a  visualization  of  the  experiment  (not  drawn  to  scale).  This  figure 
gives  a  spliced  view  of  the  test  mechanism  with  an  illustration  of  the  two  glass  layers. 
The  spores  wells  are  directly  beneath  the  sections  of  carbon  black  which  are  1500  /irn 
apart.  The  laser  heats  the  carbon  black  which  then  emits  heat  to  the  spores. 


Glass  Slide  (1mm) 

Etched  Coverslip 
(0.15mm) 


Figure  2.2:  Visualization  of  Experiment 


The  spores  were  heated  against  the  carbon  black  wells  via  a  Nd:YAG  laser  which 
was  operated  at  different  powers  in  order  to  vary  the  temperatures  experienced  by  the 
spores.  The  laser  beam  was  rastered  across  each  column  and  row  of  the  test  apparatus. 
The  heating  times  were  controlled  by  adjusting  the  raster  rate  of  the  sample  relative 
to  the  laser  beam.  Immediately  following  the  heat  treatment,  the  B.a.  spores  were 
counted  using  a  Zeiss  microscope  and  allowed  to  incubate  and  grow.  “Growth  rates 
are  less  a  function  of  the  organism  and  more  that  of  the  sample’s  original  preparation 
and  the  environmental  conditions  of  that  preparation;  temperature,  moisture  content, 
pH”  [5].  Ideal  growth  times,  based  on  the  organisms  used,  fluctuated  between  three 
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and  five  hours.  Therefore,  the  spores  were  allocated  a  growth  period  of  three  hours 
and  were  again  viewed  and  counted  via  the  Zeiss  microscope. 

Now  armed  with  an  appreciation  of  the  general  experiment,  it  is  necessary  to 
understand  which  parameters  of  the  experiment  are  to  be  modeled.  As  one  of  the 
research  objectives  is  to  evaluate  the  rate  of  thermal  diffusion  propagated  throughout 
the  B.a.  spores,  the  numerical  model  will  only  take  into  consideration  the  thermal 
heating  of  one  spore  well.  The  model  will  illustrate  the  thermal  heating  throughout 
the  coverslip  and  slide  and  will  take  into  account  the  rate  of  heat  diffusion  from  the 
center  of  one  spore  well  to  approximately  halfway  in  between  this  spore  well  and  the 
adjacent  spore  well.  The  boxed  region  shown  in  Figure  2.3  illustrates  the  model  area 
of  concern  (not  drawn  to  scale). 
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Figure  2.3:  Area  of  Experiment  to  be  Modeled 


2.2  Model  With  Respect  to  Experiment 

2.2.1  Heat  Equation. 

This  research  is  concerned  with  the  rate  at  which  the  B.a.  spores  directly  be¬ 
neath  the  carbon  black  are  heated  and  the  influence  this  heating  has  on  adjacent 
spores.  In  other  words,  does  the  adjacent  spore  realize  any  significant  thermal  energy 
from  the  previously  heated  carbon  black  well?  To  answer  these  questions,  we  consider 
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the  thermal  diffusion  equation  as  given  by: 


pC 


dT 

~di 


V  •  KVT  +  S 


(2.1) 


where 

T  =  temperature 
t  =  diffusion  time 
C  =  specihc  heat 
p  =  density 

K  =  thermal  conductivity 
S  =  a  distributed  heat  source. 

In  relation  to  the  experiment,  our  source  has  a  power  density  proportional  to  the 
Nd:YAG  laser  intensity  and  its  units  are  W/cm2.  See  Table  2.1  for  the  properties  of 
soda  lime  glass  used  in  the  experiment  [42], 


Table  2.1:  Properties  of  Soda  Lime  Glass 


Property 

Value 

Units 

C 

P 

K 

K  = 

pc 

.879 

2.44 

0.94xl(T2 

0.438xl0-2 

J/g°c 

g/ cm3 
W/cm°C 
cm2 /sec 

Because  the  spore  wells  and  the  carbon  black  are  uniformly  spaced  and  each 
is  circularly  symmetric,  we  can  consider  the  model  area  to  be  a  three  dimensional 
cylindrical  volume.  As  such,  we  will  be  evaluating  the  heat  equation  in  cylindrical 
coordinates.  Further,  we  assume  the  material  is  homogeneous,  isotropic,  and  has 
no  angular  dependence.  Thus,  the  model  is  reduced  to  axis  symmetric  and  two  di¬ 
mensional  with  respect  to  the  radial  direction  and  throughout  the  thickness  of  the 


2-4 


apparatus.  Therefore,  Equation  (2.1)  becomes: 


dT  Id,  dT.  d2T  S 
~  KrWr(r  ad +  ks37  +  ^ 


(2.2) 


where  r  is  the  radial  distance  and  z  is  the  distance  through  the  thickness  of  the  test 
apparatus. 


2.2.2  Carbon  Black. 

Lamp  black  has  an  emissivity  of  e  =  0.98.  As  such,  we  will  assume  carbon  black 
to  be  a  perfect  black  body  absorber.  A  black  body  absorbs  all  light  that  falls  onto 
it  and  emits  as  much  thermal  energy  as  it  absorbs  [10].  No  light  passes  through  a 
black  body  and  none  is  reflected.  These  properties  make  the  region  of  carbon  black  an 
ideal  source  of  thermal  energy.  In  the  experiment,  the  carbon  black  region  is  the  only 
source  of  the  thermal  energy  that  heats  the  region  of  the  spores  through  conduction. 
Refer  to  Figure  2.4(d)  to  view  the  carbon  black  immediately  after  heat  treatment. 

2.2.3  Source  Term. 

The  heat  source  is  produced  by  the  Nd:YAG  laser  beam.  The  beam  has  a 
circular  Gaussian  distribution  with  a  half-power  beam-width  of  250  /jrn  in  diameter. 
Therefore,  if  the  beam  is  centered  on  an  absorber  well,  the  heat  source  is  nearly 
uniform  across  the  carbon  black  well.  The  laser  provides  the  radiant  energy  which  the 
carbon  black  absorbs.  Due  to  experimental  limitations,  the  carbon  black  wells  were 
not  populated  equally  with  carbon.  This  reduction  in  carbon  black  diminishes  the 
efficiency  of  a  given  carbon  black  well.  Refer  to  Figure  2.4  for  microscopic  views  of  the 
carbon  black  wells.  Only  the  wells  that  were  filled  85%  (see  Figure  2.4(b))  to  100%  (see 
Figure  2.4(c))  with  carbon  black  were  used  as  viable  wells  for  the  experiment.  Since 
the  experiment  only  took  into  account  those  carbon  black  wells  that  were  populated 
by  at  least  85%,  and  because  the  wells  were  only  75  /jrn  wide,  the  laser  beam  was 
sufficiently  large  enough  to  assume  a  uniform  illumination.  This  leads  to  uniform 
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(a)  Carbon  Black  well  that  was  barely  filled  (b)  Carbon  black  well  filled  by  at  least  85%. 
and  thus  not  considered  a  viable  well  for  the 
experiment . 
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(c)  Perfectly  filled  carbon  black  well. 

(d)  Completely  filled  carbon  black  well  viewed 
after  heat  treatment. 

Figure  2.4:  Carbon  black  wells  as  viewed  underneath  the  Zeiss  microscope  [5]. 
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source  in  the  carbon  black  region  proportional  to  the  laser  power  and  the  volume  of 
the  deposited  carbon.  Further,  the  source  is  assumed  to  either  be  on  or  off  and  only 
exists  in  the  carbon  black  region.  The  reason  for  this  is  that  the  glass  is  assumed 
to  be  a  perfectly  non-absorbing  medium,  hence  a  totally  ineffective  radiator.  This 
is  a  good  assumption  for  our  temperatures  of  interest.  For  example,  the  difference 
between  hot  and  cool  glass  at  600  °C  cannot  be  seen  but  it  can  be  seen  at  1000  °C. 
For  this  reason  we  model  the  source  as 


S(t) 


kP,  laser  on 
0,  laser  off 


where  P  is  the  laser  power  in  Watts  and  k  is  the  proportionality  constant  (this  is 
determined  during  laboratory  calibration). 


2.2.4  Boundary  Conditions. 

The  boundary  conditions  need  to  imitate  the  laboratory  experiment  as  closely 
as  possible.  As  shown  by  Figure  2.3,  the  left-hand  boundary  of  the  model  area  is  the 
center  of  the  carbon  black  and  spore  wells.  By  the  axisymmetric  assumption,  there 
will  be  no  heat  flux  across  this  boundary  as  heat  must  flow  outward  from  the  wells. 
Therefore,  we  will  impose  the  condition 

=  0.  (2.3) 

r= 0 

The  test  apparatus  is  held  together  by  clamps  at  the  corners  of  the  slides  and 
suspended  in  air  to  allow  rastering  of  the  laser.  Therefore,  at  the  top  and  the  bottom 
of  the  slide,  the  boundary  condition  allows  for  convective  heat  transfer  through  air.  In 
other  words,  the  “flux  across  the  surface  is  proportional  to  the  temperature  difference 
between  the  surface  and  the  surrounding  medium”  [9].  We  will  use  Newton’s  law  of 
cooling  to  describe  this  heat  transfer  from  the  surface  of  the  slides  to  the  surrounding 
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air.  This  leads  to  the  equations 


and 


where 


-  =  KT  -  T.)  |«=,op 

OT 

—  ^s)  lz=2bottom 


Ts  =  ambient  temperature 
h  =  surface  convection  constant. 


(2.4) 

(2.5) 


For  our  model  we  can  assume  Ts  =  0  and  compute  the  change  in  temperature  with 
respect  to  ambient  temperature. 

It  is  convenient  to  define  the  dimensionless  parameter  Bi  called  the  Biot  number. 
The  Bi  number  is  a  measure  of  efficiency  of  the  heat  transport  across  the  boundary 
with  regard  to  internal  heat  conductance.  It  is  defined  by 


For  small  values  of  the  Biot  number,  the  test  apparatus  is  weakly  coupled  to  the 
ambient  temperature  outside  the  glass.  For  large  values  of  the  Biot  number,  there  is 
a  strong  coupling  to  the  ambient  temperature  and  the  surface  tends  to  approach  this 
outside  temperature. 


2.2.5  Characteristic  Scaling. 

Certain  characteristic  scaling  was  introduced  to  convert  the  variables  into  di¬ 
mensionless  variables.  This  allows  our  calculations  to  be  unit  free  or  independent  of 
particular  units  chosen.  In  other  words,  the  physical  law  (i.e.  the  heat  equation)  re¬ 
mains  valid  regardless  of  the  units  chosen  to  express  the  variables  [6] .  “The  solutions 
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of  problems  in  conduction  of  heat  can  always  be  expressed  in  terms  of  a  number  of 
dimensionless  quantities”  [9] .  The  dimensionless  variables  are  defined  in  the  following 
manner: 


r  =  Iv, 


z  =  K, 


12t 


where 


T  =  uT0, 


and 


T0  = 


SI 2 

~Y 


l  =  length 

v  =  dimensionless  length 
(  =  dimensionless  thickness 
r  =  dimensionless  time 
u  =  dimensionless  temperature 
Tq  =  initial  temperature. 


Before  proceeding  with  the  numerical  evaluation  that  utilizes  these  characteristic 
scales,  the  assumptions  formulated  for  the  thermal  model  will  be  discussed. 


2.2.6  Assumptions. 

Due  to  constraints  involved  with  modeling  a  laboratory  experiment,  certain 
assumptions  were  made  during  development  of  the  numerical  model.  The  assumptions 
are  summarized  below. 

•  The  source  term  is  uniform  and  is  considered  to  be  either  on  or  off. 

•  The  carbon  black  absorbs  all  optical  energy  and  converts  it  to  thermal  energy. 

•  The  glass  absorbs/radiates  none  of  the  thermal  energy. 

•  The  carbon  black  region  has  the  same  thermal  properties  as  the  glass  slides. 

•  Certain  characteristic  or  dimensionless  scales  are  used  in  the  model. 
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In  addition,  one  parameter  of  the  laboratory  experiment  described  above  was 
altered  after  completion  of  this  research’s  thermal  model.  The  current  width  of  the 
coverslip  wells  is  50  fin i  which  was  increased  from  3  fm 1.  The  depth  of  the  wells 
was  also  increased  from  1  //in  to  a  varying  3  to  5  //m.  This  expanded  well  size  and 
subsequent  population  of  the  wells  with  spores  allowed  air  to  be  trapped  within  the 
wells.  The  addition  of  air  in  between  the  spores  and  the  soda  lime  glass  is  significant 
because  only  the  glass  properties  were  considered  in  the  original  thermal  diffusion 
model,  and  heat  dissipates  through  glass  and  air  at  a  different  rate.  But  diffusivity, 
k,  for  soda  lime  glass  is  0.438  x  10-2  cm2/sec  while  k  =  0.2216  cm2/sec  for  air  so  the 
ratio  of  diffusivity  (^air/^giass)  is  50.559.  This  implies  that  the  rate  of  heat  propagation 
through  the  test  apparatus  with  the  addition  of  the  air  pocket,  is  much  smaller  than 
the  rate  of  the  original  model.  Since  we  are  primarily  concerned  with  the  thermal  effect 
on  the  adjacent  spore  well  and  the  original  model  overestimates  this  measurement, 
this  parameter  was  not  altered.  These  assumptions  and  the  characteristic  scales  were 
used  to  develop  and  implement  the  dimensionless  numerically  evaluated  model. 


2.2.7  Numerical  Methods. 

By  utilizing  characteristic  scaling,  the  axisymmetric,  two-dimensional  thermal 
diffusion  system  of  equations  (Equations  (2.2),  (2.3),  (2.4),  and  (2.5))  becomes 


du  d2u  1  du  d2u 
dt  dr 2  r  dr  +  dz2 


du  . 

—  (0,z,t)  =  0 

du .  .  . 

-Q^{r,zt0p  ,t)  =  -Biu 

du  . 

v  ’  ^'bottom?  tj  Bill 
u(r,  z,  0)  =  0 


(2.6) 

(2.7) 

(2.8) 
(2.9) 

(2.10) 


where  t  is  time,  r  is  the  radius,  z  is  the  distance  through  the  thickness  of  the  test 
apparatus,  and  s  is  the  dimensionless  source. 
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This  system  of  equations  is  numerically  evaluated  via  a  second-order  Crank- 
Nicolson  finite-difference  scheme  with  Strang  splitting.  A  second-order  method  im¬ 
plies  the  truncation  error  propagated  by  the  Crank-Nicolson  scheme  is  on  the  order  of 
(At)2  where  At  is  the  time  step  chosen  for  the  model.  More  importantly,  the  Crank- 
Nicolson  method  has  many  stability  advantages.  It  is  an  implicit  method  and  it  is 
absolutely  stable  which  implies  that  for  a  stable  differential  equation  and  step  size, 
any  perturbation  created  in  one  time  step  is  not  amplified  in  subsequent  time  steps. 

Before  numerically  discretizing  this  equation,  we  implement  Strang  splitting  by 
splitting  the  A  operator  into  its  dimensional  components  and  apply  half  of  the  source 
to  each  dimension.  Thus,  Equation  (2.6)  is  converted  into 

du*  d2u  1  du  x 

dt  dr2  ^  r  dr  2<S 

du  d2u*  1 

dt  =  ~d^  + 

Strang  splitting  is  applied  in  order  to  maintain  the  second  order  error  caused  by 
splitting  the  operator  in  time.  It  is  implemented  by  solving  Equation  (2.11)  for  half 
of  a  time  step,  then  solving  Equation  (2.12)  for  a  full  time  step,  and  finally  solving 
Equation  (2.11)  for  another  half  time  step  [25]. 

Let  U  be  the  numerical  approximation  to  u,  specifically 

Uj:k  =  U  (r  j  i  tn) 


(2.11) 

(2.12) 


where 


rj  =  (j  —  l)Ar,  Zk  =  kAz,  and  tn  =  nAt 

and  A r  is  the  spatial  step  in  the  radial  direction  and  Az  is  the  spatial  step  in  the  z 
direction  (thickness)  while  n,  j,  and  k  are  the  number  of  iterations  such  that  1  <  j  < 
N  and  1  <  k  <  N.  Through  use  of  the  Crank-Nicolson  discretization  with  respect  to 
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space  and  time,  Equation  (2.11)  becomes 


-  Vb  mtb.  -  2U"i' '  +  vpt*)  +  (t^M  -  +  E?-u)] 


At 


4r'j  Ar 


[(c 


n+ 1 

1+1,  k 


At 


U-+lk)  +  (£/?+u  -  £?_,,*)]  +  (2.13) 


Similarly,  Equation  (2.12)  leads  to 


Tjn+l 

Uj,k 


At 

2Az2 


~  2 UV+1  +  U%\) 


At 


+  (U-k+1-2U-k  +  U-k_1)]+-s 


(2.14) 


Our  model  has  a  uniform  mesh  or  equal  grid  spacing  in  both  the  radial  and  thickness 
directions,  that  is,  A r  =  A z. 

Along  the  left-hand  boundary  of  the  model  area  (where  r  =  0),  the  cylindrical 
heat  equation  becomes  invalid  due  to  the  term  in  Equation  (2.11).  This  issue 
is  resolved  by  treating  t  as  a  constant  and  expanding  U  around  r  —  0  via  Taylor 
series  expansion.  By  means  of  this  method,  a  centered-difference  approximation  can 
be  obtained  and  Equation  (2.7)  becomes 


T-rn+i  _  jjn  _ 
Ul ,fc  Ul,k  ~ 


2A  t 
(A  r)2 


mr  -  os') + wb  -  pw] 


(2.15) 


where  j  —  1  is  the  left-hand  boundary  [31].  Equation  (2.15)  is  applied  within  the 
model  at  the  center  of  the  left  side  of  the  domain  (where  r  =  0).  Near  the  right  side  of 
the  model  area  we  will  be  measuring  the  effect  of  the  thermal  heating  as  it  propagates 
toward  the  adjacent  spore  well.  Dirichlet  boundary  conditions  are  implemented  at  the 
right-hand  domain  far  enough  from  our  sampling  that  it  does  not  impact  the  solution. 
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The  boundary  conditions  along  the  top  and  bottom  of  the  slides  are  also  found 
by  means  of  a  difference  approximation.  Equations  (2.8)  and  (2.9)  give 

U-.tr1  ~  Ujj,  =  Btl.At  (Uft1  -  l®')  (2.16) 

(®l  -  u"i  =  Bibo«o (£/”«  -  (/”,)  (2.17) 

where 

Bhop  2  +  BiAz2  d  ^bottom  2  —  BiAz2 ' 

The  top  boundary  of  the  slides  is  defined  to  be  k  =  N  while  k  —  1  is  the  bottom 
boundary.  Equations  (2.16)  and  (2.17)  are  applied  within  the  model  at  the  top  and 
bottom  of  the  domain.  The  thermal  model  consisting  of  Equations  (2.13),  (2.14), 
(2.15),  (2.16),  and  (2.17)  was  implemented  and  numerically  evaluated  in  MATLAB 
in  order  to  determine  the  effects  of  thermal  diffusion  throughout  the  spores  and  the 
test  apparatus.  For  details  of  the  MATLAB  code,  refer  to  Appendix  A. 

2.3  Results 

As  discussed  in  Section  2.2.1,  our  model  is  concerned  with  two  main  measure¬ 
ments.  The  first  is  a  comparison  of  the  temperature  that  the  spores  are  experiencing 
as  opposed  to  the  temperature  that  the  carbon  black  is  encountering.  These  results 
are  illustrated  in  Figure  2.5.  The  top  line  is  the  dimensionless  temperature  calculated 
at  the  center  of  the  source  (the  carbon  black  well).  The  bottom  line  represents  the 
dimensionless  temperature  that  the  spores  directly  beneath  the  center  of  the  source 
are  subjected  to.  The  spore  temperature  is  calculated  4  /mi  below  the  carbon  black 
which  corresponds  to  the  average  depth  of  the  spore  wells.  The  x-axis  represents 
time  and  clearly  shown  in  the  plot  is  a  0.12  second  application  of  the  thermal  source. 
The  0.12  second  exposure  reflects  the  fastest  raster  rate  utilized  within  the  laboratory 
experiment.  The  y-axis  is  dimensionless  temperature,  or  to  be  more  exact,  it  is  the 
percentage  temperature  change  from  ambient  temperature.  By  comparing  the  respec¬ 
tive  dimensionless  temperatures  at  0.12  seconds,  the  spores  were  found  to  experience 
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approximately  93.1%  of  the  heat  that  the  carbon  black  emitted.  Also  of  note  is  the 
time  it  takes  for  the  temperatures  of  the  spores  and  carbon  black  to  ramp  up  to  a 
steady  state  value.  Both  the  spores  and  the  carbon  black  reach  approximately  92% 
of  their  total  dimensionless  temperature  within  0.014  seconds.  They  reach  85%  at 
approximately  6.85  x  10~3  seconds.  These  measurements  were  calculated  with  a  heat 
transfer  coefficient  of  h  =  2  x  10-3W/cm2oC  which  led  to  Bi  =  .0213  [43]. 


Figure  2.5:  Source  and  Spore  Temperature  Profile  for  0.12  Second  Exposure 

The  second  parameter  of  interest  within  the  model  is  the  temperature  effect  felt 
by  the  adjacent  spore.  Figure  2.6  illustrates  the  dimensionless  temperature  profile  at 
the  point  halfway  in  between  spore  wells  (refer  to  Figure  2.3).  This  characterization 
was  also  calculated  for  a  0.12  second  application  of  the  source.  It  takes  approximately 
0.05  seconds  for  any  heat  to  dissipate  from  the  carbon  black  well  to  the  point  of 
interest.  The  max  temperature  is  realized  at  approximately  0.18  seconds  which  is 
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0.06  seconds  after  the  source  is  turned  off.  Most  significant  is  that  the  temperature 
experienced  at  this  point  is  on  the  order  of  10-4  less  compared  to  the  maximum 
temperature  experienced  at  the  center  of  the  carbon  black.  Therefore,  since  this 
profile  is  taken  halfway  in  between  adjacent  source  wells,  the  effect  of  a  carbon  black 
well  on  the  adjacent  spore  is  negligible. 


Figure  2.6:  Temperature  Profile  Halfway  Between  Adjacent  Carbon  Black  Wells 


2.4  Conclusion 

The  thermal  model  developed  in  this  chapter  produces  two  main  results.  The 
change  in  temperature  experienced  by  the  spores  is  approximately  93.1%  of  the  tem¬ 
perature  change  seen  by  the  carbon  black.  This  measurement  allows  a  more  accurate 
development  of  a  probability  of  kill  model  for  the  thermal  heating  method  designed  in 
the  laboratory  experiment.  More  important  is  the  depiction  of  the  temperature  pro¬ 
file  halfway  between  adjacent  carbon  black  wells.  The  negligible  effect  of  one  carbon 
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black  well  on  the  adjacent  spore  is  also  necessary  to  develop  an  accurate  probability 
of  kill  model.  In  addition,  this  calculation  allows  modification  of  future  experiments 
in  order  to  ease  the  challenges  that  the  1500  /jrn  separation  of  carbon  black  wells 
created  within  the  laboratory.  The  characterization  of  thermal  heating  of  a  spore 
raises  the  question  of  exactly  how  the  spore  is  being  damaged.  The  next  two  chapters 
describe  our  proposal  of  the  events  occurring  within  a  spore  that  causes  the  spore  to 
be  damaged  to  such  an  extent  that  the  spore  is  effectively  dead. 
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III.  Diffusion  Model 


3. 1  Introduction 

Now  that  our  method  of  thermal  exposure  is  understood,  we  will  examine  the 
effects  of  the  thermal  damage  on  the  spore  itself.  By  applying  heat  to  a  spore, 
we  energize  adsorbed,  absorbed,  and  chemically  bound  water  molecules.  Adsorbed 
water  molecules  are  bound  to  the  outside  surface  of  a  spore  and  they  become  mobile 
when  heated  to  temperatures  between  approximately  100  to  120  °C.  Absorbed  and 
chemically  bound  water  molecules  are  found  within  the  spore.  Absorbed  molecules 
become  available  at  temperatures  similar  to  adsorbed  water  while  chemically  bound 
water  molecules  become  available  only  at  temperatures  above  approximately  300  °C. 
These  measurements  are  consistent  with  thermo-gravimetric  analysis  studies  which 
confirm  significant  weight  losses  of  Bacillus  spores  when  heated  to  temperatures  of 
100  °C  and  300  °C  [22],  Some  chemically  bound  water  molecules  are  released  through 
Maillard  reactions.  Maillard  reactions  are  chemical  reactions  between  sugars  and 
proteins  that  are  initiated  by  heat.  Maillard  reactions  occur  within  different  phases 
but  two  of  these  phases  cause  water  molecules  to  break  away  from  proteins  or  amino 
acids. 

We  propose  energizing  of  water  molecules  allows  them  to  become  free  water 
molecules  which  causes  the  molecules  to  be  more  active  in  hydrolysis  reactions  with 
DNA  and  proteins,  and  permits  mobility  or  diffusion  of  the  water  within  the  spore. 
In  addition,  diffusion  of  water  between  the  spore  and  the  environment  surrounding 
it  also  occurs.  This  interaction  between  the  spore’s  outside  environment  is  related 
to  the  water  concentration  of  the  environment.  For  example,  during  dry  heating  the 
water  concentration  of  the  air  surrounding  the  spore  would  depend  on  the  ambient 
humidity.  During  a  wet  heating  process,  the  water  concentration  of  the  surrounding 
liquid  would  effect  the  diffusion  of  water  in  and  out  of  the  spore.  The  release  and 
mobility  of  water  leads  to  hydrolysis  reactions  within  the  spore  which  will  be  discussed 
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in  further  detail  within  Chapter  IV.  In  order  to  accomplish  this,  we  must  first  have 
an  appreciation  of  the  structure  of  a  spore. 

3.2  Structure  of  Bacillus  Spores 

Bacillus  spores  are  constructed  of  concentric  layers  all  of  which  have  different 
roles  in  regards  to  the  spore’s  heat  resistance.  The  layers  of  a  Bacillus  spore  are 
illustrated  in  Figure  3.1  (layers  not  drawn  to  scale).  We  will  discuss  each  layer  in 
detail. 


CELL 
WALL 


INNER 

MEMBRANE 


OUTER 

MEMBRANE 


Figure  3.1:  Structure  of  Spore  [40] 


The  outermost  layer  is  the  exosporium  and  as  already  stated,  the  exosporium 
is  not  present  on  all  species  of  the  Bacillus  genus.  Most  importantly,  the  exosporium 
is  not  believed  to  provide  extra  protection  or  heat  resistance  [21].  The  exosporium  is 
mostly  made  up  of  proteins. 

Inside  the  exosporium  lies  the  spore  coat  which  is  composed  of  several  layers. 
It  is  also  made  up  of  proteins.  It  has  been  shown  that  heat  resistance  does  not 
depend  on  the  spore  coat  [16].  Immediately  inside  the  spore  coat  lies  the  outer 
membrane.  Although  the  exact  function  of  the  outer  membrane  is  not  understood, 
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the  outer  membrane  “has  no  notable  effect  of  spore  resistance  to  heat”  [40].  In 
addition,  it  is  not  thought  to  be  a  significant  permeability  barrier.  Westphal  and 
others  found  that  an  increase  in  relative  humidity  caused  a  spore  to  swell  due  to 
water  intake  [46].  Therefore,  we  will  assume  that  the  exosporium,  the  spore  coat,  and 
the  outer  membrane  are  all  permeable  with  respect  to  water  molecules. 

The  next  layer  of  a  spore  is  the  cortex,  which  is  made  up  of  peptidoglycan.  The 
cortex  reduces  the  water  content  of  the  spore  core  by  allowing  small  molecules,  like 
water,  to  pass  through.  Similarly,  the  cortex  keeps  DNA,  which  lies  in  the  core,  dry. 
The  germ  cell  wall  which  is  also  made  up  of  peptidoglycan  becomes  the  cell  wall  of 
the  outgrowing  spore  once  germination  occurs.  The  germ  cell  wall  provide  no  heat 
resistance  for  the  spore  [20,40]. 

The  inner  membrane  separates  the  core  from  the  cortex  and  it  is  believed  that 
the  inner  membrane  provides  a  strong  permeability  barrier  to  protect  the  core  from 
any  damage.  However,  Xu,  Labuza,  and  Diez-Gonzalez  heated  B.a.  spores  within 
sodium  hypochlorite  and  found  that  this  caused  the  leakage  of  dipicolinic  acid  (DPA). 
Since  DPA  is  contained  within  the  core,  this  study  proved  one  method  of  increasing  the 
inner  membrane’s  permeability  [47].  Another  study  by  Coleman,  Chen,  Li,  Cowan, 
and  Setlow  demonstrated  a  release  of  DPA  during  a  moist-heat  treatment  which 
“suggested  that  the  permeability  barriers  that  block  movement  of  small  molecules  in 
and  out  of  the  spore  core  had  been  breached”  [11].  Therefore,  this  research  assumes 
a  certain  permeability  of  the  inner  membrane  as  the  spore  membranes  are  damaged. 

The  innermost  layer  of  a  spore  is  the  core.  The  core  contains  the  DNA,  ribo¬ 
somes,  most  spore  enzymes,  and  certain  proteins.  Depending  on  the  species,  water 
makes  up  27  to  55%  of  the  core  wet  weight.  In  addition,  the  amount  of  free  water 
is  extremely  low  which  indicates  that  molecular  mobility  within  the  core  is  almost 
nonexistent.  The  core  water  content  is  the  major  contributor  to  a  spore’s  resistance 
to  wet  heat.  A  spore  is  more  resistant  to  wet  heat  with  a  low  core  water  content  [19]. 
Also  significant  is  the  fact  that  core  water  content  is  dependent  upon  sporulation,  or 
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reproduction,  temperatures  and  spore  preparation  temperatures.  Core  water  content 
has  an  inverse  relationship  with  preparation  temperatures,  the  higher  the  prepara¬ 
tion  temperature,  the  lower  the  core  water  content  [28].  Due  to  an  understanding  of 
a  spore’s  structure  and  the  varying  roles  in  relation  to  heat  resistance,  we  can  now 
develop  a  model  to  illustrate  mobility  of  water  within  a  spore. 

3.3  Reaction  Kinetics 

Before  proceeding  with  the  mobility  model,  it  is  necessary  to  discuss  chemical 
reactions  that  will  occur  within  the  spore.  These  chemical  reactions,  which  are  known 
as  reaction  kinetics,  are  governed  by  rates  of  reactions.  The  main  factors  influencing 
the  reaction  rates  are  the  concentration  and  physical  state  of  the  reactants,  the  tem¬ 
poral  behavior  of  the  reactants,  and  the  presence  of  any  catalysts  in  the  reaction  [32], 
For  example,  in  Section  4.3  we  will  model  the  rate  of  DNA  damage  by: 

d  tDNA]  =  -h  [DNA]  -  k-2  [H20]  [DNA] 


where 

[DNA]  ^information  content  of  DNA 

k\  =rate  coefficient  associated  with  [DNA]  breakdown  during  pyrolitic 
thermal  damage 

h>2  =rate  coefficient  associated  with  [DNA]  breakdown  during  hydrolysis 
[H2O]  =average  water  concentration. 

The  first  term  represents  pyrolysis  reactions  which  damage  DNA  at  high  temperatures 
while  the  second  term  represents  hydrolysis  reactions  which  damage  [DNA]  at  lower 
temperatures  if  water  is  present.  The  units  of  [DNA]  and  [H2O]  are  given  by  molality. 
Molality  is  denoted  by  molal  and  defined  as  the  number  of  moles  of  solute  per  kilogram 
of  solvent.  In  the  case  of  [DNA]  it  represents  the  concentration  of  base  pair  encoding 
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information.  The  concept  of  reaction  kinetics  will  be  used  throughout  the  water 
mobility  and  kill  probability  models. 

3.4  Water  Mobility  Model 

The  water  mobility  model  will  be  based  on  a  production/diffusion  equation. 
Molecular  diffusion  is  the  movement  of  molecules  from  a  region  of  higher  concentration 
to  that  of  a  region  with  a  lower  concentration.  This  transport  of  molecules  is  founded 
on  random  molecular  movement.  As  we  are  modeling  the  diffusion  throughout  a 
spore,  we  have  a  three  dimensional  volume  to  consider.  Diffusion  is  mathematically 
modeled  via  the  multi-dimensional  Fick’s  law  which  predicts  how  diffusion  affects  the 
concentration  of  water  over  time  [27]: 

^  =  f{w)  +  V-DVw  (3.1) 

where 

w  =  concentration  of  mobile  water  (rnolal) 
t  =  time  (sec) 

f(w)  =  production  of  water  (molal/sec) 

D  =  diffusion  coefficient  (cm2/sec). 

The  diffusion  coefficient  is  a  measure  of  the  mobility  of  water  within  the  spore.  Note 
that  when  f(w)  is  negative,  water  is  being  consumed.  Similarly,  when  f(w)  is  positive, 
water  is  produced. 

3-4-1  Initial  Distribution. 

Due  to  the  properties  of  a  spore’s  structure,  there  are  only  two  main  spore  layers 
to  consider  when  discussing  the  initial  distribution  of  water  throughout  a  spore.  These 
layers  are  the  core,  which  contains  DNA,  and  the  cortex.  As  previously  stated,  the  core 
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contains  an  extremely  low  amount  of  liquid  water  molecules,  or  absorbed  water.  Thus 
for  the  purpose  of  our  model,  we  will  consider  the  concentration  of  absorbed  water  in 
the  core  to  be  zero.  But  taking  into  account  that  water  makes  up  a  certain  percentage 
of  the  core  wet  weight,  we  must  assume  chemically  bound  water  exists  within  the  core. 
The  cortex  allows  mobility  of  water  molecules  so  it  contains  absorbed  water  molecules 
along  with  its  chemically  bound  water  molecules.  Refer  to  Figure  3.2  for  a  graphical 
illustration  of  a  notional  distribution  of  water  throughout  the  spore.  As  we  will  be 
scaling  the  radius  of  a  spore  to  one,  the  radius  of  the  core,  rc,  is  approximately  75%  of 
the  radius  of  the  entire  spore,  designated  by  rQ  [34,46].  Note  that  we  are  assuming  an 
equal  amount  of  chemically  bound  water  throughout  the  core  and  cortex  while  only 
the  cortex  contains  absorbed  water.  The  distribution  of  water  outside  the  spore  will 
depend  on  ambient  water  concentration,  which  in  turn,  depends  on  whether  a  wet  or 
dry  heating  process  is  utilized. 


Figure  3.2:  Initial  Distribution  of  Water  Throughout  Radius  of  Spore 
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3-4-2  Release  of  Water. 

Now  that  we  have  determined  the  initial  distribution  of  available  water  through¬ 
out  the  spore,  we  need  to  calculate  how  this  absorbed  and  chemically  bound  water 
is  released.  In  other  words,  how  does  the  water  become  mobile?  We  propose  that 
by  applying  heat  to  the  spore,  the  absorbed  and  chemically  bound  water  molecules 
become  energized.  Once  energized,  the  water  becomes  mobile  in  a  free  molecular  state 
similar  to  a  gaseous  state. 


[H20],  +  Heat  — *•  [U20}g 
[H20]fe  +  Heat  — ►  [H20]fl 


where 


[H20]^  =  initial  concentration  of  absorbed  water  molecules  (molal) 

[H20]6  =  initial  concentration  of  chemically  bound  water  molecules  (molal) 
[H20]  =  initial  concentration  of  mobile  water  molecules  (molal). 


These  two  sources  of  mobile  water  are  energized  at  different  temperatures.  Ab¬ 
sorbed  water  is  released  at  approximately  100  °C  while  bound  water  requires  a  much 
higher  temperature  of  approximately  300  °C.  We  choose  to  model  their  respective 
release  functions  as  a  Gaussian  density  function  centered  about  their  respective  mean 
release  temperatures  Tf  and  T2  with  associated  standard  deviations  <7i  and  cr2,  i.e. 


9i(T) 


i  =  1, 2. 


These  functions  are  illustrated  in  Figure  3.3  using  Ti=100  °C,  cr1=20,  T2=300  °C, 
and  cr2=50.  We  observe  that  absorbed  water  is  released  in  a  neighborhood  of  100°C 
until  all  the  absorbed  water  is  active.  Then  bound  water  begins  its  release  to  a  peak 
of  300  °C.  These  measurements  were  chosen  to  approximate  the  thermo-gravimetric 


3-7 


analysis  results  found  by  Holwitt,  Kiel,  Alls,  Morales,  and  Gifford  [22],  Their  study 
found  that  approximately  10%  of  the  initial  weight  of  a  Bacillus  thuringiensis  spore 
was  lost  by  100  °C  and  40%  was  lost  by  approximately  300  °C. 


Figure  3.3:  Distribution  of  Water  Released  at  100  and  300°C 


The  total  mobile  water  released  at  any  given  temperature,  T*,  is  determined  by 
integrating  the  weighted  sum  of  these  Gaussian  functions  i.e. 


[H20]9  =  f  [H2O]j0i(t)  +  [H20]b^2(r)dr 

J — oo 

=  [H20](/I(T-)  +  [H20]i,/2(T-) 


(3.2) 


where 


1  +  erf 


i  —  1,2  (3.3) 


and  erf  is  the  error  function  [1],  We  will  assume  that  given  a  particular  T* ,  this 
temperature  is  reached  instantaneously  and  therefore  the  initial  amount  of  released 
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water  in  the  spore  is  dependent  on  T* .  This  is  consistent  with  Chapter  II  results  that 
illustrated  the  immediate  increase  in  temperature  once  the  heat  source  was  applied 
(refer  to  Figure  2.5). 

Figure  3.4  shows  the  results  of  this  integration  using  the  values  from  Figure  3.3, 
[H20]£=1,  and  [H20]b=3.  This  figure  illustrates  the  effect  on  the  available  mobile 
water  as  temperatures  increase  through  the  release  temperature  for  absorbed  water 
to  the  release  temperature  for  bound  water. 


Figure  3.4:  Total  Water  Released  as  Temperature  Increases 


This  function  serves  as  our  initial  distribution  of  mobile  water  in  the  spore  for 
the  water  mobility  model.  That  is, 


w(r,0)  =  [H20]s(r) 


(3.4) 


where  r  is  the  radius. 


3-9 


3-4-3  Analytical  Analysis. 

3. 4-3.1  Production  Term. 

Water  is  consumed  through  reaction  with  hydrocarbon  compounds  within  the 
spore  such  as  DNA  and  proteins,  thus  the  production  term  in  Equation  (3.1)  is  reduced 
to  a  simple  first  order  reaction  term: 

f(w)  =  —ktw  (3.5) 

where  kt  is  the  total  reaction  rate  constant.  It  is  effectively  the  averaged  rate  for  all 
possible  hydrolysis  reactions  in  the  spore.  The  value  kt  follows  the  Arrhenius  equation 
which  models  the  rate  of  chemical  reactions  [27]: 

kt{T)  =  Ae~^r 


where 


A  =  encounter  frequency 
E  =  activation  energy 
R  =  gas  constant. 

3. 4-3. 2  Boundary  Conditions. 

For  the  purpose  of  this  analysis,  we  assume  a  spore  is  a  symmetric  sphere.  We 
model  the  diffusion  of  water  radially  from  the  spore’s  core  into  the  cortex  and  the 
surrounding  environment.  At  the  boundary  where  the  radius  is  equal  to  zero,  the 
solution  must  be  finite,  that  is, 


w(0,  t) |  <  oo. 


(3.6) 
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Where  the  radius  is  equal  to  the  total  spore  radius  (r  =  ra),  the  boundary  con¬ 
dition  follows  a  similar  argument  to  the  Newton’s  law  of  cooling  boundary  condition 
found  in  Section  2.2.4  of  the  thermal  model.  In  other  words,  the  flow  of  water  across 
the  outside  boundary  is  proportional  to  the  water  content  in  the  spore  at  the  outer 
boundary  and  the  ambient  humidity.  The  boundary  condition  is  modeled  by, 

-  n  ■  D\7w  (r G,  t)  =  h{w  (r0,  t)  -  [H20]a}  (3.7) 


where 


n  =  outward  normal  vector 
[H20\a  =  water  concentration  outside  spore  (molal) 
h  =  surface  flow  velocity. 

Note  that  we  are  assuming  the  external  [H20]a  vapor  is  at  the  spore  temperature 
T*.  Also  for  spherical  geometry,  n  becomes  the  outward  radial  direction.  Therefore, 
Equation  (3.7)  becomes 

-  D ^  (r  0,  t)  =  h{w  (r0,  t)  -  [H20]a}  .  (3.8) 

The  combination  of  Equations  (3.1),  (3.2),  (3.3),  (3.4),  (3.5),  (3.6),  and  (3.8)  consti¬ 
tute  the  water  mobility  model. 

3. 4-3. 3  Analytical  Solution. 

We  assume  the  spore  is  spherically  symmetric,  that  is,  it  has  no  angular  depen¬ 
dence.  This  allows  the  model  to  become  one  dimensional  with  respect  to  the  radius 
of  the  spore  and  time  ( w(r,t )).  Further,  we  shall  assume  the  diffusivity,  D,  remains 
constant  across  the  radius  of  the  spore.  Thus,  the  diffusion  term  in  Equation  (3.1) 
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becomes 


V  •  D\7w 


D 


d2w  2dw\ 
dr 2  r  dr  ) 


(3.9) 


At  this  point  we  introduce  dimensionless  variables  by  choosing  characteristic 
scales  for  length,  time,  and  water  concentration.  To  this  end,  we  scale  all  lengths 
by  the  radius  of  the  spore  ra,  time  by  the  diffusion  time  tc  =  r2/D ,  and  water 
concentration  by  the  maximum  concentration  of  initial  released  water  [H20]M,  i.e. 


[HjO]^  max  [H20]  (r,  T*). 

0  <r<r0  y 

We  will  assume  in  this  model,  that  [H20]a  and  [H20]6  are  constants  and  define  di¬ 
mensionless  variables  x,  r,  and  v  by 


r  =  rQx,  t  —  tcr ,  w  (r ,  t)  =  [H20]M  v  (x,t)  .  (3.10) 


Substituting  these  variables  into  Equations  (3.1),  (3.2),  (3.6),  and  (3.8)  leads  to 


dv 

=  -0v  +  -^-(x2—) 
x2  dx  \  dx) 

0<X<1,  T>0 

(3.11) 

dr  .  . 

&(0'r) 

=  0 

T  >  0 

(3.12) 

dv 

dx ^  ,T^ 

=  B(v(1,t)  -7) 

T  >  0 

(3.13) 

v  (x,  0) 

=  a1{x)f1{T*)  +  a2{x)f2{T*) 

0  <  x  <  1 

(3.14) 

where 


with 


and 


[H20L 

[H20]m 


a\(x) 


[H2O], 

[H20]m 


(■ rQx ) 


and 


a2(x) 


[H20]& 

[H20]m 


(■ rax ). 


These  parameters  have  important  physical  interpretations.  First,  we  observe 
that  7  is  the  ratio  of  the  water  concentration  at  the  boundary  of  the  spore  to  that 
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within  the  spore.  For  wet  heating  we  would  expect  7  to  be  greater  than  one  while  for 
dry  heating  this  parameter  would  be  less  than  one.  Next,  we  examine  B  and  /3  as  the 
ratio  of  time  scales.  There  are  three  time  scales  which  naturally  present  themselves 
in  this  water  diffusion  model.  They  are  the  diffusion  time,  tjj,  the  hydrolysis  reaction 
time,  tR,  and  the  surface  convection  time,  tg-  They  can  be  defined  as  follows: 

rn  1  1  ro 

tD  =  77,  tR  =  — ,  and  ts  =  — 

D  kt  h 

recalling  that  kt  is  the  total  reaction  rate  constant.  These  times  can  be  ordered  from 
smallest  to  largest,  for  example  tR  <  to  <  ts-  This  ordering  suggests  the  priority 
at  which  these  processes  take  place.  So  for  the  above  ordering,  hydrolysis  reactions 
occur  much  faster  than  diffusion  through  the  spore  or  convection  away  from  the  spore 
surface.  Further,  by  definition  we  see 


B 


tD 

ts 


1 


and  (3  =  —  >  1 
t-R 


with  this  ordering.  Because  diffusivity,  D ,  surface  convection,  h,  and  the  reaction 
coefficient,  kt ,  are  temperature  dependent,  this  time  scale  ordering  may  change  with 
temperature  and  influence  the  mobile  water  concentration  found  in  hydrolysis  reac¬ 
tions. 

Using  separation  of  variables,  the  homogeneous  problem  given  by  Equation 
(3.11)  and  the  homogeneous  boundary  conditions  becomes 


h'{r)  +  (P  +  n)h{r)  =  0 

(3.15) 

[. x2g'(x )]7  +  /j,x2g(x)  =  0 

(3.16) 

</(o)  =  0 

(3.17) 

-9'{l)  =  Bg{l) 

(3.18) 

where  p  is  a  separation  parameter.  The  eigenvalues  and  eigenfunctions  are  found  by 
solving  Equations  (3.16),  (3.17),  and  (3.18)  with  respect  to  the  homogeneous  bound- 
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ary  conditions.  The  associated  eigenfunctions  are  of  the  form 


gn(x)  =  K , 


sin(snx) 


n  e  N 


where 


AT  = 


I  2(1  —  B) 
sin2(s„)  —  B 


is  found  by  normalizing  gn  and  sn  satisfies  the  transcendental  equation 


tan(sn)  = 


1  —  B' 


Notice  these  eigenfunctions  could  also  be  represented  in  terms  of  spherical  Bessel 


functions  as  [1] 


9n{x)  AnSnJo  (Sn*^)  • 


These  eigenfunctions  are  orthonormal  with  respect  to  the  weight  x2,  i.e. 


(fln?  9m)  /  %  Qn(x)Qrn(x) d.X 


Now  we  seek  a  solution  to  the  homogeneous  system  of  the  form 


v(x,  t)  =  7  +  ^2  hn{r)gn{x). 


(3.19) 


Introducing  Equation  (3.19)  back  into  Equation  (3.11)  gives 


h'n(r)  +  l unhn(r)  =  -/3y  (1,  gn) 


(3.20) 


where 


wn  —  +  s2 


3-14 


Equation  (3.20)  along  with  the  initial  conditions  (Equation  (3.14))  produces  the  so¬ 
lution  to  the  first-order  differential  equation: 

hn{r)  =  Fn(T*)e-“nT  -  [p  +  s2ne~WnT]  (3.21) 

where  Fn(T*)  =  fi(T*)  (a1}gn)  +  f-iiT*)  (a2,gn). 

Inserting  Equation  (3.21)  into  Equation  (3.19)  and  reducing  terms  gives 

oo  oo  2 

v{x, T)  =  Fn(T*)e-“nTgn(x)  +  7  (1,  gn)  —  [l  -  e_aJ"r]  gn(x).  (3.22) 

1  1  ^71 

n=  1  7i=l 

The  Erst  series  represents  the  influence  of  the  initial  water  distribution  while  the 
second  series  represents  the  influence  of  the  wet  or  dry  heating  (denoted  by  7).  Recall 
from  Equation  (3.10)  that 


w(r,t)  =  [H20 ]mv(x,t)  . 


(3.23) 


DNA  is  only  present  in  the  spore  core.  Thus  all  the  hydrolysis  reactions  with 
DNA  will  occur  in  the  core  so  we  are  only  concerned  with  the  diffusion  of  water  within 
the  spore  core.  The  average  total  mobile  water  concentration  in  the  core  is  given  by 


[H20]c  (t)  =  f  w(r,t)r2dr.  (3.24) 

rc  Jo 

Equation  (3.24)  with  the  inclusion  of  the  dimensionless  variables  and  Equation  (3.22), 
reduces  to 


[H20]c  (ter)  =  4  [H20]m  \Y,K(T‘)Gn(xe)e-“- 

y  n=l 

00  2 
+7^Gn(l )G„.(a)-  [l  —  e-w”T] 


71=1 


(3.25) 
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where  xc  is  defined  as  rc/r0  and 


Gn{x)=  r^gniO  de  (3.26) 

Jo 

Notice  if  x  —  1  then  Equation  (3.26)  gives 

Gn{  1)  =  [  x2gn{x)dx  =  (1  ,gn) . 

Jo 

Equation  (3.25)  contains  important  limiting  behaviors.  For  short  time  periods 
(small  values  of  t),  the  second  term  approaches  zero  implying  that  the  first  term,  or 
the  influence  of  local  water,  dictates  the  release,  diffusion,  and  reaction  of  water  in 
the  core.  Only  during  longer  time  periods  does  the  effect  of  the  water  concentration 
outside  the  spore  become  relevant.  For  more  details  of  the  analytical  solution,  see 
Appendix  B. 


3. 4- 3. 4  Steady  State  Solution. 

The  steady  state  solution  is  obtained  from  Equation  (3.11)  where  is  assumed 
to  vanish.  Along  with  the  boundary  conditions  (Equations  (3.12)  and  (3.13)),  the 
steady  state  problem  becomes 


_d_ 

dx 


—  f3x2vs  =  0 


dvs 

dx 

dus 

dx 


(o)  =  o 

(1)  =  B(  v,(l)-y) 


where  s  indicates  steady  state.  This  boundary  value  problem  admits  a  solution 


A 


vs(x)  =  —  sinh(v//fa) 
x 
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where 


A  = 


7  B 


cos1i(a//3)  —  (1  —  B)  sinh(v/^) 

We  are  assuming  that  the  environment  cannot  be  perfectly  dry  thus  7  ^  0.  Now 
using  Equations  (3.23)  and  (3.24),  the  steady  state  core  water  concentration  becomes 


[H20]«=-r[H20]Mt.s 


X 


(3.27) 


with 


X  an 


7 B  T  \[flxc  cos1i(a//Tcc)  —  sinh(A/^a:c) 

/ 3  _v//3cosh(v//3)  -  (1  -  B)  sinh(V^) 


(3.28) 


The  steady  state  solution  represents  the  limiting  behavior  of  the  time  dependent 
core  water  concentration  given  in  Equation  (3.25).  In  fact,  Equation  (3.25)  can  be 
written  as 


[H20]c  (tcr)  =  [H20]sc  +  -  [H20]m  £  Cne 


n=  1 


where 


a  = 


Fn(T*)  —  7Gn(l)  — 


Ur, 


Steady  state  is  reached  at  approximately  half  the  diffusion  time  ( ts  ~  ^to)- 


3.5  Water  Mobility  Model  Results 

3. 5. 1  Assumptions. 

During  development  of  the  water  diffusion  model  which  examines  water  release, 
diffusion,  and  reaction  within  the  spore  core,  the  following  assumptions  were  made: 

•  The  geometry  of  the  spore  and  its  properties  are  spherically  symmetric. 

•  The  initial  distribution  of  water  depends  on  adsorbed,  absorbed,  and  chemically 
bound  water  molecules  outside  and  throughout  the  spore. 

•  These  water  molecules  become  gaseous  water  molecules  once  released  by  heating 
the  spores  to  a  certain  temperatures. 
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•  The  desired  temperature,  T*\  which  is  reached  instantaneously. 

•  Diffusion  can  only  take  place  with  gaseous  water  molecules. 

•  The  gaseous  water  molecules  are  consumed  in  chemical  reactions  with  other 
molecules  such  as  DNA  and  proteins  (hydrolysis). 

In  addition  to  these  assumptions,  the  parameters  found  in  Table  3.1  were  used 
to  solve  the  production/diffusion  model  [8,14,30,33,48]. 


Table  3.1:  Water  Mobility  Model  Parameters 


Parameter 

Value 

Units 

R 

D 

h 

A 

E 

[H20]m 

1.986x  10-3 
2.0xl0~9 
1.023x10s 
1.6xl03 

10 

25.944 

kcal/moPK 
cm2  /sec 
cm/sec 
sec-1 
kcal  /  mol 
molal 

3.5.2  Aiialytical  Results. 

Utilizing  the  assumptions  and  parameters  described  above,  the  model  was  solved 
analytically  to  evaluate  the  influence  that  time,  temperature,  and  heating  process  (wet 
versus  dry)  had  on  the  spore  core  water  content.  Figure  3.5  was  analyzed  for  7  =  1  at 
a  temperature  of  110  °C.  This  7  value  indicates  that  the  water  concentration  outside 
the  spore  was  equivalent  to  the  initial  released  gaseous  water  concentration  inside 
the  spore.  The  time  behavior  of  the  core  water  concentration  is  indicated  by  the 
solid  line  while  the  steady  state  result  (refer  to  Equations  (3.27)  and  (3.28))  is  given 
by  the  dashed  line.  As  previously  discussed,  over  time  the  core  water  level  will  rise 
as  water  is  released  and  will  equilibrate  according  to  its  relationship  with  the  water 
concentration  outside  the  spore.  Using  the  values  in  Table  3.1,  the  diffusion  time,  to, 
is  1.25  seconds  and  steady  state  is  reached  in  approximately  0.6  seconds. 
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Figure  3.5:  Core  Water  Concentration  for  7=1  at  110  °C 


Figure  3.6:  Core  Water  Concentration  for  7=1  at  320  °C 
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Figure  3.6  maintains  the  same  value  for  7  (7  =  1)  but  is  analyzed  at  an  initial 
release  temperature  of  320  °C.  When  compared  to  Figure  3.5,  the  long  term  core 
water  level  is  essentially  equivalent  for  both  temperatures  but  the  difference  in  initial 
water  concentration  is  caused  by  the  release  of  chemically  bound  water  due  to  the 
higher  temperature. 

Now  we  will  examine  the  evolution  of  the  core  water  level  when  a  spore  is 
thermally  treated  via  a  wet  heating  process.  Figure  3.7  shows  the  change  in  core 
water  concentration  at  320  °C  for  7  =  10.  This  implies  that  the  water  concentration 
outside  the  spore  is  much  higher  than  inside  the  spore  (wet  heating).  Note  that 
the  amount  of  water  in  the  core  increases  rapidly  to  a  much  higher  level  as  water  is 
diffusing  into  the  cortex  and  then  into  the  core  from  outside  the  spore. 


Figure  3.7:  Core  Water  Concentration  for  7=10  at  320  °C 


Figure  3.8  demonstrates  the  evolution  of  the  core  water  concentration  for  7  = 
0.1,  or  a  dry  heating  process,  at  320  °C.  This  example  clearly  illustrates  the  limiting 
time  behaviors.  Initially  the  core  water  level  rises  since  all  the  released  absorbed 
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and  chemically  bound  water  in  the  cortex  immediately  flows  outward.  Keep  in  mind 
that  the  core  does  not  contain  any  absorbed  water  so  the  cortex  water  concentration 
initially  seeks  to  equilibrate  by  diffusing  into  the  core.  But  as  time  progresses,  the  ratio 
of  ambient  humidity  to  spore  water  concentration  becomes  the  dominating  behavior 
and  water  flows  out  of  the  spore.  In  addition,  the  core  water  content  level  is  much 
lower  than  the  previous  examples. 


Figure  3.8:  Core  Water  Concentration  for  7=0.1  at  320  °C 


As  a  comparison,  Figure  3.9  graphs  the  core  water  concentration  for  the  same 
value  of  7  but  at  an  initial  water  release  temperature  of  110  °C.  Since  the  lower 
temperature  gives  a  constant  value  of  released  absorbed  water  within  the  spore,  there 
is  no  initial  water  diffusion  from  the  cortex  into  the  core,  just  an  equilibrium  with  the 
surrounding  environment . 

The  release,  diffusion,  and  reactions  of  water  within  a  spore  cause  damage  to  a 
spore’s  DNA  and  proteins.  The  results  of  the  water  mobility  model  and  subsequent 
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Figure  3.9:  Core  Water  Concentration  for  7=0.1  at  110  °C 


damage  due  to  hydrolysis  reactions  will  now  be  utilized  to  characterize  a  spore’s 
probability  of  kill. 
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IV.  Probability  of  Kill  Model 


4-1  Introduction 

Chapter  III  illustrated  the  release  and  mobility  of  water  throughout  a  spore. 
This  diffusion  of  water  allows  hydrolysis  reactions  to  occur  within  the  spore.  As 
explained  in  the  introduction,  we  contend  that  hydrolysis  reactions  leading  to  de- 
ploymerization  of  the  spore’s  DNA  and  proteins  will  decimate  these  long  chain  molecules 
to  such  an  extent  that  they  can  no  longer  be  repaired.  We  further  assume  that  damage 
to  the  spore  is  based  on  a  certain  initial  protein  fitness  of  the  spore  population.  Once 
this  fitness  level  is  degraded  beyond  a  critical  level,  the  spore  can  no  longer  repair 
damaged  DNA  and  thus  cannot  germinate,  or  produce  outgrowth.  These  assertions 
allow  us  to  create  a  probability  of  kill  model  based  on  water  mobility,  hydrolysis, 
DNA  information  content,  and  fitness  levels  of  the  spore’s  proteins. 

4-2  Hydrolysis 

Once  water  molecules  become  mobile  and  diffuse  throughout  a  spore,  chemical 
reactions  occur.  Hydrolysis  is  a  reaction  of  a  molecule  with  water.  The  acid  or 
base  environment  for  hydrolysis  is  important  in  the  following  way.  Hydrolysis  can  be 
catalyzed  by  H+  or  by  OH-,  so  often  rates  of  reaction  can  be  represented  as: 

Rate  =  (kn)  [Biomolecule] 

where  kB  is  the  total  reaction  rate  given  by 

kH  =  kA  [H+]  +  kB  [OH-]  +  kN 

where 

kN  =  k*N  [H20] 
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with 


[•]  =  activity  of  ions 
kA  =  acidic  rate 
ks  =  basic  rate 
kjsr  =  neutral  rate. 

The  measurement  of  pH  in  a  solution  is  related  to  its  acidity  or  basicity.  Bacteria 
can  survive  an  environment  with  pH  measurements  between  five  and  eight.  A  solution 
with  a  pH  below  seven  is  acidic  and  a  pH  above  seven  is  basic.  The  rate  of  reaction 
of  hydrogen  molecules  versus  pH  level  is  characterized  by 

pH  =  -  log  [H+] 

and  the  rate  of  reaction  with  pH  values  close  to  seven  is  shown  Figure  4.1  [12]. 


(a)  Log  rate  versus  pH  of  H+.  (b)  Log  rate  versus  pH  of  OH  . 

Figure  4.1:  Log  Rate  of  Reaction  versus  pH  Level 


The  relationship  between  hydroxide  molecules  is  related  to  the  pH  level  by 

—  log  Kw  —  pH  +  pOH 
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where 


Kw  =  dissociation  constant  ~  10  14  at  25  °C. 

Therefore  a  neutral  solution  that  is  neither  acidic  or  basic  has  a  concentration  of 
hydrogen  and  hydroxide  ions  equal  to  10-' moles  dm-3.  Note  that  the  rates  of  reaction 
increase  as  the  activity  of  negatively  or  positively  charged  ions  increases  [12].  More 
importantly,  “microbial  cells  require  their  internal  pH  to  remain  constant  in  order 
to  maintain  essential  cell  functions”  [41].  We  assume  that  kA  and  ks  are  constant 
corresponding  to  a  neutral  acid/base  environment  and  therefore  k n  depends  only  on 
the  rate  of  water  activity. 

As  previously  stated,  water  can  react  with  biopolymers  such  as  DNA  and  pro¬ 
teins  resulting  in  depolymerization  or  removal  of  a  side  group  of  the  polymer.  De- 
ploymerization  is  breaking  a  large  strand  of  polymers  into  two  smaller  ones.  Water 
can  also  react  with  DNA  causing  deamination  and  depurination  damage.  Deamina¬ 
tion  takes  place  when  an  amino  group  from  the  DNA’s  base  pair  is  severed.  This 
prevents  the  DNA  from  replicating  properly.  “Thus  a  change  in  the  DNA  base  code 
would  both  propagate  and  replicate  as  the  essential  cell  proteins  would  be  improperly 
produced  or  not  produced  at  all  because  of  this  DNA  damage”  [20]. 

DNA  is  most  susceptible  to  depurination  which  occurs  when  DNA  undergoes 
spontaneous  decomposition  in  living  cells.  During  depurination,  a  DNA  base  pair 
is  removed  from  the  DNA  chain.  When  DNA  is  missing  one  of  its  base  pairs,  the 
DNA  strand  can  not  match  up  with  its  other  half  of  the  helix  chain.  Information 
is  encoded  in  the  DNA  polymer  as  a  series  of  base  pairs.  DNA  contains  the  genetic 
information  enabling  organisms,  including  bacteria,  to  replicate.  The  information  is 
encoded  by  the  location  of  bases  on  the  DNA  backbone.  Each  of  four  bases  along 
the  DNA  backbone  represents  a  bit  of  information.  Three  bits  encode  for  a  letter 
in  the  DNA  language,  called  a  codon.  A  specific  codon  codes  for  a  particular  amino 
acid  in  its  place  along  the  length  of  a  protein  polymer  during  its  synthesis.  The 
information  is  preserved  and  replicated  by  the  Watson-Crick  base  pairing  in  which 
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complementary  bases  recognize  each  other.  One  incorrect  amino  acid  can  dramatically 
alter  the  function  of  the  protein,  particularly  shape  sensitive  proteins  like  enzyme 
catalysts.  Three  bits  of  ordered  information  allow  64  possible  encoding  combinations. 
The  capability  of  the  three-bit  letter  is  more  than  sufficient  to  encode  for  20  amino 
acids  that  make  up  proteins.  Thus  there  exists  redundancy  in  the  code.  As  a  result, 
removing  or  altering  the  third  base  pair  in  the  sequence  of  three  bits  in  a  codon  may 
or  may  not  destroy  the  information  content;  removing  the  first  base  causes  problems 
more  frequently.  Some  codon  letters  are  inactive,  so  destroying  that  information 
causes  no  degradation  in  function  of  the  DNA  or  the  proteins  for  which  it  codes. 
Hence,  damage  to  a  particular  DNA  base  can  have  a  range  of  effects  depending  on 
the  information  content  of  the  codon  or  the  information  content  of  the  bit  within 
the  codon.  We  limit  ourselves  to  consider  damage  which  reduces  information  content 
in  critical  or  essential  codons,  which  contain  information  necessary  to  replicate  the 
spore  during  germination.  Even  so,  there  is  a  range  of  information  degradation  that 
can  occur  for  each  base  that  is  damaged  or  removed  [13].  We  can  account  for  this 
variability  in  our  model  by  including  a  distribution  of  damage  influences  for  each  base 
that  reacts  by  hydrolytic  depurination  or  deamination. 

Water  reactions  with  proteins  cause  denaturation  damage.  Denaturation  is  a 
process  in  which  proteins  lose  their  structure  by  application  of  some  external  stress. 
Native  structures  of  proteins  may  be  altered,  and  their  biological  activity  changed 
or  destroyed  by  treatment  that  does  not  disrupt  the  primary  structure.  Following 
denaturation,  some  proteins  will  return  to  their  native  structures  under  proper  condi¬ 
tions.  Extreme  conditions,  such  as  strong  heating,  usually  cause  irreversible  change. 
Since  the  structure  of  the  protein  determines  its  function,  the  protein  can  no  longer 
perform  its  function  once  it  has  been  denatured  [37] .  For  the  purpose  of  our  research, 
we  will  assume  that  hydrolysis  reactions  are  equally  likely  to  occur  at  all  protein 
sites.  Further,  we  will  assume  that  the  critical  threshold  damage  to  DNA  caused  by 
depurination  is  equally  likely  for  all  base  pair  sites.  The  damage  caused  by  hydrolysis 
reactions  will  be  considered  in  the  probability  of  kill  model. 
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4-3  Probability  of  Kill  Model 


There  exist  many  models  describing  the  inactivation  of  Bacillus  spores.  These 
models,  through  use  of  various  heating  methods,  relate  time  and  temperature  to  the 
inactivation  of  spores.  Refer  to  Figure  4.2  for  a  graphical  representation  of  some  of 
these  models  [17,44]  where  heating  time  is  a  log  linear  function.  The  data  gathered  by 
AFIT  researchers,  Goetz  and  Hawkins,  used  similar  heating  techniques  as  described 
in  Chapter  II  [20,21],  Battelle  used  a  series  of  furnaces  to  suspend  spores  in  heat  for 
various  periods  of  time  [2], 


Log  Heating  Time  (sec) 


■ 

B.a.  S  potts- Whitney,  et  al. 

■ 

B.a.  Ferneliusetal. 

■ 

B.a.  Goetz 

■ 

B.a.  Hawkins 
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B.t.  Hawkins 

• 

B.t.  Battell  5  decades 

+ 

B.t.  Battell  2  decades 

— ©- 

HE  Region 

o 

Incendiary  Region 

Figure  4.2:  Heat  Inactivation  of  Bacillus  Spores 
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The  Battelle  data  is  the  closest  comparison  found  in  the  literature  of  short 
duration,  high  temperatures  similar  to  the  parameters  utilized  in  this  research.  Notice 
that  the  slope  of  the  heat  inactivation  of  the  Battelle  data  is  different  than  the  slope 
of  inactivation  when  spores  were  heated  at  lower  temperatures  and  for  longer  time 
periods. 

Our  model,  while  still  taking  into  account  the  effects  of  heating  temperature  and 
duration,  will  also  consider  a  third  variable  to  explain  the  inactivation  of  spores.  We 
suggest  that  diffusion  of  water  throughout  a  spore,  and  thus  the  hydrolysis  reactions 
with  DNA  and  proteins,  damage  these  molecules  to  such  an  extent  that  the  spore  is 
effectively  dead.  In  addition,  these  hydrolysis  reactions  cause  the  difference  in  the 
heat  inactivation  slope  when  spores  are  heated  at  high  temperatures  and  for  short 
durations  (refer  to  Figure  4.2)  [7].  We  assume  that  although  hydrolysis  reactions 
are  equally  likely  to  occur  at  all  DNA  sites,  these  hydrolysis  reactions  can  occur 
multiple  times  causing  the  DNA  to  lose  information  content  per  hydrolysis  event. 
Thus  we  propose  that  DNA  information  content  is  based  on  the  effect  of  damage  due 
to  hydrolysis  reactions.  We  further  propose  that  a  population  of  heat-treated  spores 
has  a  certain  initial  ‘fitness’  distribution  based  on  viability  of  protein  and  capacity  to 
repair  DNA,  and  by  choosing  one  spore  from  the  sample  we  are  provided  a  specific 
protein  fitness  level.  We  offer  a  model  which  correlates  spore  inactivation  with  time, 
temperature,  an  initial  DNA  information  content  level  within  a  particular  spore,  and 
an  initial  protein  fitness  distribution  of  the  spore  population. 

Let  In(t)  be  the  DNA  information  content  in  a  given  spore  at  time  t.  The  rate 
of  DNA  damage  caused  by  heating  and  hydrolysis  will  be  modeled  by 

-  ^  =  h(T)ID  +  k :2(T)  [H20]c  (t)ID  (4.1) 
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where  [H20]c  (t)  is  defined  by  Equation  (3.24)  and 
t  =time 

T  ^temperature 

k\  =rate  coefficient  associated  with  [DNA]  breakdown  during  pyrolysis 
/c2  =rate  coefficient  associated  with  [DNA]  breakdown  during  hydrolysis. 

This  equation  characterizes  the  rate  of  change  in  the  DNA  information  content 
of  a  spore.  The  first  term  is  the  first  order  breakdown  process  associated  with  direct 
thermal  pyrolysis  reactions  of  DNA  and  the  second  term  captures  the  influence  of 
hydrolysis.  The  rate  coefficients,  ki(T ),  follow  the  Arrhenius  equation  [27]: 

ki(T)=Aie~ ™  (4.2) 

where  is  the  respective  encounter  frequency  while  E%  is  the  respective  activation 
energy.  By  integrating  with  respect  to  time,  Equation  (4.1)  becomes 

ID{t)  =  Ioe-lkit+k2H(t)] 

with 

m  =  I  [H20]c  (r)dr. 

Jo 

By  using  the  core  water  concentration  model  in  Equations  (3.25)  and  (3.26)  this 
becomes 

3  f  00  1 

H(t)  =  -j  [H20]m  {  £F„(r)Gn(Z')—  [1  -  e-""‘] 

i  k'Tl 

c  ^  n=  1 

oo  2 

+7  Gn(l )Gn(xc)—j  \ujnt  —  (l  —  e  u,n*)] 

„_i  A 
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At  time  t  —  0,  the  initial  DNA  information  content,  I0,  will  be  determined  by  a 
spore’s  DNA  information  content  prior  to  heat  treatment  and  subsequent  hydrolysis. 
Since  the  DNA  of  one  spore  is  nearly  identical  to  that  of  another,  the  information 
content  will  be  represented  by  a  constant  function.  The  information  content  will  be 
degraded  due  to  its  hydrolytic  susceptibility.  “The  death  of  the  individual  cell  occurs 
when  the  number  of  functional  ribosomes  has  dropped  below  a  critical  threshold  level, 
beyond  which  the  cell  cannot  recover”  [41],  The  same  can  be  said  for  DNA.  “While 
there  can  be  repair  of  DNA  damage. . .  if  too  much  damage  has  been  accumulated 
during  spore  dormancy,  this  damage  can  overwhelm  the  capacity  of  repair  systems 
and  lead”  to  spore  death  [40] .  Thus  we  let  A  denote  the  critical  threshold  such  that  if 
Id  >  Ik,  then  the  spore’s  DNA  has  enough  viable  information  that  it  can  be  repaired 
and  thus  produce  a  cell  during  germination.  See  Figure  4.3  for  an  illustration  of  a 
spore’s  DNA  information  content  prior  to  any  heat  treatment  (t  —  0). 


Content  of  DNA  Information  within  a  Spore 

Figure  4.3:  Spore  DNA  Information  Content  at  t  —  0 

Let  Id(t)  be  defined  as  the  DNA  information  content  at  time  t  =  r.  As  time 
progress,  the  probability  that  critical  damage  is  done  to  a  spore’s  DNA  information 
content  increases.  We  will  assume  this  increase  in  critical  damage  follows  a  Gaus¬ 
sian  density  function  with  a  mean  of  Id  and  standard  deviation  designated  as  cjd- 
The  density  function  represents  the  hydrolytic  susceptibility  of  DNA  information  con- 
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tent.  Over  time,  &d  increases  as  the  likelihood  of  damage  to  an  essential  codon  also 
increases.  Figure  4.4  represents  the  degradation  in  DNA  information  content  due 
to  hydrolysis  reactions  as  time  evolves  from  t  —  0  to  t  —  T\  then  t  —  T2,  and  Id 
approaches  its  critical  kill  threshold,  I k- 
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Content  of  DNA  Information  within  a  Spore 

Figure  4.4:  Spore  DNA  Information  Content  as  it  Evolves  over  Time 

The  probability  density  function  for  critical  codon  damage  will  not  be  considered 
in  this  research.  Instead  we  assume  any  damage  is  critical  codon  damage.  This  leads 
to  the  conclusion  that  the  spore  is  killed  if  Id  <  Ik  and  will  survive  if  Id  >  h- 
However,  in  addition  to  DNA  damage,  there  is  damage  to  proteins  during  heating 
and  hydrolysis.  Since  proteins  are  necessary  to  repair  DNA,  a  damaged  spore  will 
only  be  able  to  germinate  with  the  protein  that  survives  the  heating  process.  So 
DNA  must  be  repaired  before  a  viable  cell  is  produced  during  germination  and  there 
must  be  a  suitable  amount  of  repair  protein  available  to  accomplish  this.  Let  Jp(t)  be 
the  protein  information  content  for  the  sample  population  at  time  t.  The  evolution 
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of  protein  content  is  also  given  by 


~^f  =  HT)Jp  +  h(T)  [H20]c  ( t)JP 


(4.4) 


where 

k\  =rate  coefficient  associated  with  protein  breakdown  during  pyrolysis 
k-2  =rate  coefficient  associated  with  protein  breakdown  during  hydrolysis 

and  the  rate  coefficients  are  also  modeled  by  the  Arrhenius  equation  (as  in  Equation 
4.2). 

At  time  t  —  0,  the  initial  protein  information  content,  J0,  will  be  determined 
by  the  spore  chosen  from  a  population  with  a  certain  protein  fitness  distribution 
based  on  various  preparation  parameters.  As  previously  discussed,  spore  preparation 
temperature  effects  the  core  water  content  which  in  turn  influences  a  spore’s  heat 
resistance.  In  addition,  “spores  prepared  at  different  temperatures  have  a  number 
of  significant  differences  in  properties”  besides  heat  resistance  [19,28].  Therefore 
due  to  the  multitude  of  preparation  parameters  that  affect  the  properties  of  a  spore, 
our  model  will  consider  that  a  spore  population  will  have  a  specified  protein  fitness 
distribution.  We  will  model  a  spore’s  initial  protein  content  via  a  Gaussian  density 
function.  At  time  t  —  0,  Jo,  is  randomly  drawn  from  this  function  which  has  a  mean 
of  Jo  and  a  standard  deviation  denoted  by  ap.  See  Figure  4.5  for  an  illustration  of 
the  protein  fitness  density  function  of  a  sample  population  of  spores  prior  to  any  heat 
treatment  (t  —  0).  After  integration,  Equation  (4.4)  becomes 

JP(t )  =  J0e-N+^{t)] 
with  Hit)  given  by  Equation  (4.3). 

Proteins  are  needed  to  repair  damaged  DNA.  For  this  reason,  we  are  only  con¬ 
cerned  with  repair  proteins  in  the  core  because  of  the  location  of  a  spore’s  DNA. 
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Figure  4.5:  Sample  Spore  Population’s  Protein  Fitness  at  t  —  0 

But  proteins  can  also  be  degraded  during  heating  and  hydrolysis  reactions  and  “this 
damage  can  overwhelm  the  capacity  of  repair  systems  and  lead  to”  spore  death  [40]. 
Let  Jk  be  the  critical  kill  threshold  such  that  when  Jp  >  J the  spore  contains  the 
necessary  repair  protein  information  content  in  order  to  repair  its  DNA. 

Since  Jo  is  a  random  variable  drawn  from  the  density  function,  Jp(t)  is  also 
random.  However,  the  standard  deviation  of  Jp(t),  o> ,  is  assumed  to  remain  constant 
while  the  mean,  J p,  evolves  with  time  since 

JP(t )  =  joe-[fcl*+fc2^h)]. 

Let  Jp(r)  be  defined  as  the  random  variable  taken  from  the  density  function  which 
progresses  to  time  t  —  r.  Over  time,  Jp  will  approach  its  critical  kill  threshold,  Jk- 
Refer  to  Figure  4.6  for  a  graphical  representation. 

In  order  to  determine  the  probability  of  kill,  P*,,  for  the  spore,  we  simply  need 
to  integrate  underneath  the  section  of  the  Gaussian  density  function  that  has  passed 
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Figure  4.6:  Sample  Population’s  Protein  Fitness  as  it  Evolves  over  Time 


the  critical  kill  threshold: 


Pk  (  Jt  Jk) 


1 


\f^2/KGp 


Following  integration  Equation  (4.5)  becomes 


(4.5) 


P»(V<A)  =  |[l+erf(QF)] 


(4.6) 


where 


Qp  — 


Jk  J  t 


(4.7) 


\/2(7  p 

This  probability  of  kill  analysis  will  be  evaluated  to  examine  the  influence  of  wet 
versus  dry  heating  on  a  spore. 
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4-4  Results 

The  parameters  used  within  the  probability  of  kill  model  are  found  in  Table 
4.1  [15,26,29,30,48], 


Table  4.1:  Parameters  for  Probability  of  Kill  Model 


Parameter 

Value 

Units 

R 

1.986  x  10~3 

kcal/mol0K 

Ax 

2.0  x  1012 

sec-1 

£1 

65.683 

kcal  /  mol 

A2 

2.564  x  1010 

sec-1 

e2 

29 

kcal  /  mol 

A1 

4.828  x  1015 

sec-1 

Ej 

38.3 

kcal  /  mol 

a2 

1.6  x  103 

sec-1 

e2 

10 

kcal  /  mol 

Figure  4.7  illustrates  the  probability  of  protein  survival  with  an  initial  release 
temperature,  T*,  of  300  °C.  The  dashed  line  represents  wet  heating  or  7  =  10.  The 
upper  line  reflects  a  dry  heating  process  where  7  =  0.1.  Note  that  it  takes  twice 
as  long  (0.1  seconds  versus  0.2  seconds)  for  the  dry  heating  process  to  achieve  a 
90%  reduction  in  protein  population  survival  (see  the  values  marked  by  a  *).  This 
supports  the  contention  that  an  increased  amount  of  hydrolysis  reactions  will  cause 
faster  protein  damage.  In  fact,  Figure  4.7  shows  that  at  0.1  seconds  wet  heating 
produced  a  10%  survival  probability  while  the  same  exposure  time  with  dry  heating 
produced  a  73%  survival  probability  (shown  by  the  o).  Further,  we  see  the  survival 
probability  curves  see  the  asymptote  to  a  constant  value  (approximately  0.25  seconds). 
Our  results  imply  that  the  entire  population  can  never  be  killed.  This  is  a  consequence 
of  modeling  the  protein  fitness  with  a  Gaussian  density  function  which  approaches 
zero  only  at  positive  and  negative  infinity.  Because  JT  is  bounded  below  by  zero,  Qp 
in  Equation  (4.7)  is  bounded  above  by  Jk/\/2<yp-  This  leads  to  a  minimum  survival 
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probability  (using  Equation  (4.6))  of 


1  —  P  _  i 

-L  -1  b  min 


-  erf  (  j}' 
\\f2op 


Logically  we  know  that  there  exists  a  finite  amount  of  information  in  the  spore  pop¬ 
ulation’s  protein  fitness  and  it  should  be  possible  to  kill  the  entire  population. 


Figure  4.7:  Probability  of  Protein  Survival  at  300  °C 


The  probability  of  protein  survival  during  a  wet  heating  process  (7  =  10)  at 
varying  temperatures  is  shown  in  Figure  4.8.  As  expected,  when  heated  at  increasingly 
higher  temperatures,  equivalent  protein  damage  takes  place  within  a  shorter  duration. 
Notice  that  even  at  200  °C  the  protein  damage  asymptote  is  reached  by  approximately 
0.6  seconds.  As  shown  in  Section  3.5.2,  this  time  duration  will  only  reflect  the  influence 
of  released  water  within  the  spore.  The  effect  of  ambient  water  concentration  is  not 
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reflected  during  these  early  time  periods.  Thus,  the  difference  in  the  probability  of 
kill  measurements  for  various  temperatures  represents  the  concentration  of  absorbed 
and  chemically  bound  water  within  the  spore.  For  example,  at  300  °C  the  release 
of  chemically  bound  water  causes  the  protein  to  be  damaged  faster  as  more  water  is 
available  for  hydrolysis  reactions. 


Figure  4.8:  Probability  of  Protein  Survival  During  Wet  Heating 

Figure  4.9  represents  a  dry  heating  process  (7  =  0.1)  for  initial  water  release 
temperatures  of  250  °C  and  300  °C.  The  disparity  between  these  two  temperatures 
reflects  the  lesser  amount  of  initial  released  water,  as  well  as  the  influence  of  ambient 
humidity.  Recall  from  Chapter  III  that  after  0.6  seconds  the  outside  water  concen¬ 
tration  dominates  the  diffusion  of  water  throughout  the  spore.  At  T*  =  250  °C,  some 
protein  survives  for  up  to  five  seconds  as  water  is  flowing  from  inside  the  spore  into 
the  surrounding  environment  and  there  is  less  water  available  for  hydrolysis  reactions 
within  the  core. 
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Figure  4.9:  Probability  of  Protein  Survival  During  Dry  Heating 


These  results  show  the  consequences  of  the  release  and  diffusion  of  water  through¬ 
out  a  Bacillus  spore.  Further,  the  production  of  water  allows  hydrolysis  reactions  to 
occur  with  the  spore’s  DNA  and  proteins  and  these  reactions  may  cause  sufficient 
damage  that  the  spore  cannot  be  repaired.  The  probability  of  a  spore  being  damaged 
to  this  extent  can  be  represented  by  a  probability  of  kill  model. 
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V.  Conclusions  and  Future  Work 


5.1  Conclusion 

This  research  contains  a  notional  characterization  of  damage  to  a  Bacillus  spore 
due  to  thermal  heating  and  hydrolysis  reactions.  These  hydrolysis  reactions  occur 
following  the  release  and  diffusion  of  water  molecules  throughout  the  spore  which  is 
effected  by  the  type  of  heat  treatment  used  (wet  versus  dry).  Hydrolysis  reactions 
damage  a  spore’s  DNA  and  proteins  and  the  extent  of  this  damage  may  result  in  the 
inability  of  the  spore’s  DNA  to  be  repaired  and  thus  not  germinate.  Further,  the 
damage  of  a  spore  can  be  represented  as  a  probability  of  kill  model. 

Many  assumptions  were  made  during  development  of  the  probability  of  kill 
model.  Throughout  the  entire  research  we  assumed  a  spore  was  spherically  symmet¬ 
ric.  In  addition,  the  rates  of  water  release,  water  diffusion,  and  water  consumption 
during  hydrolysis  reactions  were  similar  to  those  found  in  recent  studies  but  these 
rates  were  applied  within  this  research  which  has  not  been  replicated  by  laboratory 
experimentation.  The  actual  distribution  of  absorbed  and  chemically  bound  water 
molecules  within  the  spore  also  requires  further  exploration.  These  assumed  param¬ 
eters  significantly  impact  the  results  of  the  probability  of  kill  model  and  should  be 
researched  and  analyzed  in  greater  detail. 

Under  these  assumptions  the  release,  diffusion,  and  subsequent  reactions  of  wa¬ 
ter  within  the  spore’s  core  were  shown  to  be  dependent  on  time,  temperature,  and 
water  concentration  of  the  outside  environment.  The  temperature  dictated  the  initial 
amount  of  water  available  for  diffusion  and  hydrolysis  reactions  since  absorbed  water 
molecules  were  released  as  gaseous  water  molecules  at  approximately  100  °C  while 
chemically  bound  water  molecules  are  not  released  until  300  °C.  At  short  time  dura¬ 
tions  (less  than  0.6  seconds)  the  distribution  of  absorbed  and  chemically  bound  water 
within  the  spore  had  the  greatest  influence  on  the  core  water  concentration.  After  0.6 
seconds  the  surrounding  environment’s  water  concentration  played  a  significant  role 
in  the  transformation  of  the  core  water  concentration. 
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Once  the  absorbed  and  chemically  bound  water  molecules  became  available, 
their  hydrolysis  reactions  with  the  spore’s  DNA  and  proteins  caused  deploymeriza- 
tion  or  removal  of  a  side  group  of  the  polymer.  This  damage  was  modeled  within  a 
protein  probability  of  kill  model.  A  spore  populations  protein  fitness  was  degraded 
due  to  hydrolysis  reactions  and  an  assumed  critical  kill  threshold  was  utilized  to  de¬ 
termine  if  the  protein  was  damaged  to  such  an  extent  that  it  could  not  repair  the 
spore’s  damaged  DNA.  The  notional  protein  fitness  density  function  and  critical  kill 
threshold  impact  the  probability  of  kill  model  and  also  require  further  examination. 
The  probability  of  protein  survival  depended  whether  a  wet  or  dry  heating  process 
was  utilized.  A  wet  heating  process  allowed  increased  water  flow  into  the  spore  and 
a  greater  amount  of  hydrolysis  reactions  thus  protein  damage  takes  place  faster  than 
during  dry  heating  which  causes  water  to  flow  from  within  the  spore  into  the  surround¬ 
ing  environment.  Another  parameter  influencing  the  duration  required  to  damage  the 
entire  population  is  the  initial  water  release  temperature.  If  this  temperature  is  high 
enough  to  allow  release  of  chemically  bound  water  than  protein  is  damaged  more 
quickly  since  more  water  is  available  for  hydrolysis  reactions.  This  research  provides 
a  notional  characterization  of  a  spore’s  probability  of  kill  depending  on  time,  temper¬ 
ature,  and  hydrolysis  reactions  within  the  spore  but  further  analysis  is  necessary  to 
compare  this  research  to  real  data. 

5. 2  Future  Work 

During  development  of  the  production/diffusion  and  probability  of  kill  models 
contained  in  this  research,  the  following  areas  were  identified  as  requiring  further 
analysis.  Section  3.4.2  assumes  a  Gaussian  density  function  for  the  release  of  absorbed 
and  chemically  bound  water  in  a  heat-treated  spore.  This  is  a  valid  assumption  due 
to  the  Central  Limit  Theorem  which  states  that  any  population  with  a  significantly 
large  sample  size  will  have  a  sampling  distribution  that  is  approximately  Gaussian. 
But  a  Gaussian  density  function  approaches  negative  infinity  as  the  probability  of 
an  event  occurrence  approaches  zero.  It  does  not  take  into  consideration  that  if  the 
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probability  of  an  event  occurring  is  zero,  then  the  density  function  should  reflect  a 
value  of  zero.  As  spores  are  heated,  water  molecules  might  participate  in  hydrolysis 
reactions  as  soon  as  they  are  released  which  would  indicated  a  release  rate  of  zero, 
but  the  water  release  rate  can  never  be  negative.  Therefore,  further  analysis  and 
model  development  should  examine  the  effects  of  various  density  functions  for  the 
release  of  water  molecules.  In  particular,  the  Weibull  density  function  has  the  added 
benefit  that  the  density  function  is  zero  for  when  the  probability  of  event  occurrence 
is  zero.  Consequently,  if  no  water  is  released  then  the  probability  of  water  being 
released  is  zero.  Refer  to  Figure  5.1  for  a  graphical  representation  of  a  Weibull  density 
function  [45]. 


Figure  5.1:  Weibull  Density  Function 

Similarly,  in  Section  4.3,  we  assumed  a  Gaussian  density  function  for  protein 
fitness  levels  of  a  population  of  spores.  We  know  logically  there  cannot  exist  a  negative 
fitness  level  for  this  parameter.  Protein  fitness  can  be  completely  degraded  to  such 
an  extent  that  the  fitness  level  is  zero,  but  it  cannot  be  below  zero.  Thus  a  Weibull 
density  function  may  also  be  an  improved  assumption  for  these  parameter. 
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Further,  the  comparison  of  protein  probability  of  kill  to  real  data  is  difficult  as 
reaction  rate  parameters  are  not  well  known.  Another  recommendation  is  to  fit  the 
known  data  with  notional  response  behavior  motivated  by  the  model  and  examine 
the  extrapolated  model.  Guided  by  the  solution  given  in  Equation  (4.3),  we  could  use 
the  simplified  model 


H(t,  T )  =  A(T)t  +  B(T)  [1  -  e_aJt] 

to  represent  the  effect  of  hydrolysis.  The  coefficients  A(T),  B(T),  and  decay  rate  u> 
could  be  estimated  from  data.  This  rate  function  is  then  used  in  the  kill  probability 
model. 

Within  our  probability  of  kill  model,  spore  death  is  modeled  by  hydrolysis  dam¬ 
age  to  DNA  or  proteins  by  means  of  a  wet  or  dry  heating  process.  Ishihara  and  others 
found  that  spores  dehydrated  to  less  than  22%  of  their  initial  bound  water  could  not 
germinate  [23].  Depending  on  the  degree  of  degradation,  DNA  can  be  repaired  by 
proteins.  Therefore,  future  research  should  consider  a  spore  outgrowth  model  based 
on  the  combined  damage  of  DNA  and  proteins  following  heat  treatment.  If  there 
exists  sufficient  protein  that  functions  properly,  can  the  DNA  be  repaired  to  a  viable 
level  such  that  spore  replication  and  germination  occur?  And  since  this  repair  will 
involve  protein  synthesis,  will  the  repair  process  decimate  the  protein  fitness  level 
to  beyond  its  critical  kill  threshold  such  that  spore  death  occurs?  In  other  words, 
what  quantities  of  viable  DNA  and  proteins  will  allow  the  protein  to  repair  DNA 
after  experiencing  heat  damage?  This  approach  could  be  mathematically  modeled  in 
a  similar  manner  as  the  other  processes  contained  in  this  research  in  order  to  develop 
a  spore  outgrowth  model.  In  addition,  an  outgrowth  model  would  allow  improved 
comparisons  with  real  data. 

Rosenberg,  Cavalieri,  and  Ungers  propose  the  existence  of  repressor  and  anti¬ 
repressor  proteins  in  a  spore.  They  suggest  that  these  proteins  repress  or  initiate 
replication  of  DNA  respectively  leading  to  spore  death  or  germination  [35].  Rubinow 
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investigates  an  mathematical  model  of  the  accumulation  rates  of  repressor  proteins, 
anti-repressor  proteins  and  the  triggering  of  protein  synthesis  with  DNA  [36].  In 
order  to  extend  the  probability  of  kill  model  to  include  a  spore  outgrowth  model,  the 
proposed  repressor  and  anti-repressor  protein  initiatives  could  be  further  analyzed. 

Finally,  this  research  only  considers  hydrolysis  reactions  with  DNA  and  proteins. 
This  concept  could  easily  be  extended  to  other  reaction  kinetics.  The  same  analysis 
applies  if  treating  a  spore  with  chemicals  other  than  water. 
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Appendix  A.  MATLAB  Code  for  the  Thermal  Heating  Model 


The  following  appendix  contains  the  MATLAB  code  used  to  numerically  evaluate  the 
thermal  heating  model  in  Chapter  II: 

function  the siscyl indr i cal coordv7test 
clc; clear  all;clf; 

N=600 ; 

dt=. 00002;  "/Time  step 

°/0Radial  direction;  half  of  distance  between  adjacent  spores 
r=linspace(0, .8, (2*N/3)) ; 

dr=r(2)-r(l) ;  ”/„spatial  step  in  r  direction 
rmesh=length(r) ;  °/0develop  uniform  mesh 

%Z  direction;  total  width  of  glass  slide  and  coverslip 
z=l inspace (- . 2, 1 ,N) ; 

dz=z(2)-z(l) ;  '/ospatial  step  in  z  direction 
zmesh=length(z) ;  °/0develop  uniformmesh 

°/„Set  up  ut=urr+(l/rj)*ur  matrix,  where  rj=dr*j ,  Crank-Nicholson  method 
lambdar=dt/ (dr . ~2 ) ; 

q=spdiags (repmat ( [1  -2  1] , [rmesh  1] ) ,-l : 1 ,rmesh,rmesh)*lambdar ; 

UL3=dt . / ( (2*dr . ~2) .* (1 : rmesh) ) ; 

UL4=dt . / ( (2*dr . ~2) . * (0 : rmesh-2) ) ; 

UL4=  [0,UL4] ; 

q=(q-spdiags (UL3 ’ , -1 , rmesh , rmesh) ) ; 
q=(q+spdiags (UL4) , 1 , rmesh, rmesh) ) ; 

A=(speye(rmesh)+(0.5*q)) ; 

B=(speye(rmesh)-(0.5*q)) ; 
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°/0ut=urr+(l/r)*ur  boundary  condition  when  r=0 
B(1 , l)=l+2*lambdar ;  B(1 , 2)=-2*lambdar ; 

A(1 , l)=l-2*lambdar ;  A(1 , 2)=2*lambdar ; 


7»Set  up  ut=uzz  matrices,  Crank-Nicholson  method 

°/0Due  to  Strang  Splitting,  z  discretization  is  for  half  time  step 
dtz=dt/2 ; 

lambdaz=dtz/ (dz . ~2) ; 

Q=spdiags (repmat ( [1  -2  1] , [zmesh  1] ) ,-l : 1 ,zmesh,zmesh)*lambdaz; 
C=(speye(zmesh)+(0.5*Q)) ; 

D=(speye (zmesh) - (0 . 5*Q) ) ; 

Bi=.0213;  "/Biot  number  where  Bi=hl/K 

lambdaBib=(Bi*dt)/(2-Bi* (dz~2) ) ;  "/factor  for  bottom  boundary  conditions 
lambdaBit=(-Bi*dt) / (2+Bi* (dz~2) ) ;  "/factor  for  top  boundary  conditions 

"/Boundary  Conditions  for  the  uzz  matrix: 

"/Lower  BC,  based  on  Newton’s  Law  of  Cooling 
D(l,l)=l-lambdaBib;  D(l,2)=0;  C (1 , l)=l-lambdaBib;  C(l,2)=0; 

"/Upper  boundary  condition,  based  on  Newton’s  Law  of  Cooling 
D (zmesh , zmesh- 1) =0 ;  D (zmesh , zmesh) =l-lambdaBit ; 

u=zeros (zmesh, rmesh) ;  "/Initial  Conditions  (zero  heat  throughout  domain) 

"/Solve  ut=urr+(l/rj)*ur+uzz+S  via  Strang  Splitting 

for  j  =  1 : 5256 ;  "/Run  for  4000  loops  to  find  solutions  1  &  3, 

°/5256  loops  to  find  solution  2  (see  below) 

"/Set  up  source  term,  assign  constant  value  for  radius  of  source 
if  j<2628  '/Source  ’on’  for  .12  sec 

sour ce=sour cef un (zmesh , rmesh , dt ) ; 
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else 


source=zeros(zmesh,rmesh) ;  '/.Turn  source  'off’  after  .12  sec 

end 

u=zdirection(C,D,u, source) ;  °/0Solve  in  z  direction  for  half  time  step 
u=radial (A, B,u, source) ;  °/oSolve  in  r  direction  for  a  full  time  step 
u=zdirection(C,D,u, source) ;  °/„Solve  in  z  direction  for  half  time  step 
end 
end 

function  source=sourcefun(zmesh,rmesh, dt) 
constant=zeros(zmesh,rmesh) ; 

"/(Placement  of  source  with  4  Micron  depth  &  36  Micron  radius 
constant (100 , 1 : 18)=dt*20 ;  constant (101 , 1 : 18)=dt*20; 
source=constant ; 
end 

function  f irstst ep=zdirect ion (C,D ,u, Source) 

f irststep=D\ (C*u+0 . 25*Source) ; 

end 

function  secondstep=radial (A , B , u , Source) 
secondstep=(B\ (A*u’+0 . 5*Source ’ ) )  ’ ; 
end 
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Appendix  B.  Details  of  Water  Mobility  Analytical  Solution 


The  following  appendix  contains  details  of  the  analytical  solution  to  the  water  produc¬ 
tion/diffusion  model  found  in  Chapter  III.  The  original  production/diffusion  equation 
is 


=  -ktw  +  V  •  DVw. 
dt 


(B.l) 


With  one  dimension,  the  diffusion  term  in  Equation  (B.l)  becomes 


^  „„  1  d  (  9dw\ 

V  ■  DVw  =  —  —  Dr2  — 
rz  or  \  or  J 


(B.2) 


and  because  D  is  constant,  Equation  (B.2)  is  actually 


^  „  / d2w  2d w 

V  ■  DVw  =  D  [  ——  +  -  — 
or  r  or 


(B.3) 


At  this  point,  we  introduce  dimensionless  variables  by  choosing  characteristic  scales 
for  length,  time,  and  water  concentration.  We  scale  all  by  the  spore  radius,  rG; 
time  by  the  diffusion  time,  tc  =  r2/D ;  and  water  concentration  by  the  maximum 
concentration  of  the  initial  released  water,  [H20]M,  for  any  given  temperature,  T*, 
i.e.  [H20]m  =  max0<r<T.o  [H20]  (r,T*).  We  define  dimensionless  variables  x,  r,  and 
v  by 


r  =  rQx ,  t  =  tcr,  and  w  (r,  t)  =  [H20]M  v  (x,  r) .  (B.4) 


Therefore,  with  the  dimensionless  variables,  Equation  (B.3)  becomes 


V  ■  DVw 


D  [H20]m 

12  P20]M 

'  o 

W  [h20]m 

'  o 


/  d2v  1  2  dv  \ 

\  dx 2  D  xrl  dx  J 
/  d2v  2  dv  \ 

\  dx2  x  dx  ) 

1  d  f  2dv\ 
x2  dx  \  dx  J 


(B.5) 
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Inserting  Equation  (B.5)  and  other  dimensionless  variables  into  Equation  (B.l)  gives 


[HgQ \_mOv 
tc  dr 
dv 
frr 


1  d  (  2dv 

-0v  +  ^ax[xai 


where 


With  dimensionless  variables,  the  original  outer  boundary  condition, 


—  D 


dw 

dr 


(r0,  t ) 


h{w  (r0,t)  -  [H20]J 


(B.6) 


becomes 


~7„  |H2°Im  £(1't)  =  A{[H20]Ml.(l,r)  -  [H2OU 

-£(1,T)  =  S(l’(1.T)  -7) 


in  which  the  total  radius  of  the  spore  is  scaled  to  the  dimensionless  value  of  1, 


B  = 


rQh 

D  ’ 


and 


[H-.OL 

[h20]m 


In  addition,  the  initial  condition  (t  —  0), 


w(r,  0)  =  [H20]a  (r)  =  [H20],  (r)/,(T*)  +  [H20]„  (r )/2(T*) 


with  dimensionless  variables  becomes 


v(x,  0)  = 


w(rQx,  0) 


[H20]m  [H20]m 

=  «i  (i)/,(T*)  +  a2(i)/2(T*) 


[H20]^  (rp*\  ,  [H20]6  (rox}f2(T* 


r0x)f1(T*)  + 


[H20] 


M 
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where 


,  ,  [H20],  . 

“lW  =  TiTTvi  \rox) 


and 


(  S  [H20]6 

Mx)  =  ftt  ('W 


[h2o]m^'  [H20]m 

In  summary,  the  dimensionless  production/diffusion  equation  with  boundary  condi¬ 
tions  and  initial  conditions  is 

dv  1  d  (  2dv\ 

Ih  =  ~(5v+h?~dx  V  dx) 

v  (x,  0)  =  a\  (x)  fi  (T*)  +  a2(x)f2  (T*) 

To  solve  this  system  of  equations,  we  will  use  separation  of  variables.  Let 

v  ( x ,  r)  =  g  (x)  h  (r) , 

then  Equation  (B.7)  becomes 

g(x)h'(r)  =  —(3h{r)g(x)  +  |W(x)]  • 


0  <  X  <  1,  T>0 

(B.7) 

r  >  0 

(B.8) 

r  >  0 

(B.9) 

0  <  x  <  1. 

(B-10) 

Let 


ti(r)  +Ph(r)  _  yj  ([xV(x)]_ 

—  /r, 


(B.ll) 


/r(r)  3(x) 

where  ji  ia  a  separation  parameter.  Thus  Equations  (B.7),  (B.8),  and  (B.9)  become 


h\r)  +  ((3  +  /i)h(r)  =  0 

(B.12) 

^2  [^V^)] '  +  ng{x)  =  0 

(B.13) 

</(0)  =  0 

(B.14) 

-g\l)  =  Bg{l). 

(B.15) 
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To  simplify  the  solution  to  Equations  (B.13),  (B.14),  and  (B.15),  let 


g(x) 


J(x) 

x 


then 


g'(x) 


xJ'(x)  —  J(x) 


(B.16) 


or 

x2g\x )  =  xJ'(x)  —  J(x). 


Since 


x2g{x)  =  xJ(x), 


Equation  (B.13)  is  reduced  as 


[xJ'(x)  —  J(x)}'  +  / jlxJ(x )  =  0 
J\x)  +  xJ"(x)  —  J'(x )  +  gxJ(x)  =  0 
J"(x )  +  fiJix)  =  0. 


The  solution  to  this  differential  equation  is  of  the  form 


J(x)  =  Ai  cos(sx)  +  A2  sin(sx) 


where  s  =  y/Jl.  Substituting  into  Equation  (B.16), 


g\x)  = 


-Aisx  sin(sx)  +  A2sx  cos(sx)  —  A1  cos (sx)  —  A2sin(sx) 


x z 


=  —  Ais2 


sxsin(sx)  +  cos(sx) 


sx 


+  4L2s2 


sxcosfsxj  —  sm(sxj 
(sx)2 
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To  satisfy  Equation  (B.14),  linxE_»0+  g'(x)  =  0  requires  Ai  =  0.  Further,  Equation 
(B.15)  becomes 


g'(l)  +  B  g(l)  =  J\l)  -  J(l)  +  B  J(l) 

=  A2  [scos(s)  —  (1  —  B)  sin(s)] 
=  0 


For  A2  ^  0  we  require 


s  cos(s)  —  (1  —  B)  sin(s)  =  0 


or 


tan(s) 


s 

1  -B 


(B-17) 


provided  B  ^  1.  If  B  —  1  then  we  require  cos(s)  =  0  or  sn  =  (2 n  —  1)|.  This  same 
result  is  found  in  the  limit  as  B  — >  1  in  Equation  (B.1T).  Therefore,  the  eigenvalues 
and  eigenfunctions  associated  with  Equations  (B.13),  (B.14),  and  (B.15),  with  respect 
to  the  homogeneous  boundary  conditions,  are  of  the  form 


gn(x)  =  K , 


sin(sna;) 


tie  N 


x 


(B.18) 


where  the  constant  Kn  is  found  by  normalizing  gn  and  the  eigenvalues,  sn ,  are  de¬ 
termined  by  satisfying  Equation  (B.17).  Figure  B.l  represents  a  graphical  solution 
to  this  equation  for  values  of  B  equal  to  0.2  and  5.  It  shows  the  intersection  of  the 
line  y  —  s/(l  —  B)  with  the  tangent  function.  The  points  of  intersection  determine 
the  eigenvalues.  Note  that  the  eigenvalues  for  B  =  0.2  are  indicated  by  a  +  sign  on 
the  x-axis  while  a  o  sign  represents  the  eigenvalues  for  B  =  5.  Further,  we  see  as  n 
increases,  the  value  of  sn  approaches  (2 n  —  1)7t/2.  In  fact  a  good  approximation  for 
large  values  of  n  is 


sn  ~  (2 n 


1) - tan  1 

;  2 


1  -B  ' 
_(2n  —  1)|.  ' 
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Figure  B.l:  Eigenvalues  For  B  Equal  to  0.2  and  5 

Table  B.l  presents  the  first  eight  values  of  sn  for  various  values  of  B  between  0.2  and 
10.  Note  as  n  increases,  the  eigenvalue  approaches  (2 n  —  l)vr/2  and  the  sequence  con¬ 
verges  to  this  value  faster  for  large  values  of  B.  Figure  B.2  graphically  illustrates  the 
values  of  s  with  n  equal  to  1,  2,  and  3  as  B  varies  from  0  to  10.  These  eigenfunctions 
are  orthonormal  with  respect  to  the  weight  x 2,  i.e. 

(dm  9rn)  /  •£  9n  (x}gm  (ic)diC  8 mn • 

J  0 

To  prove  that  gn  and  gm  are  orthogonal  we  use  the  ordinary  differential  equation 
(Equation  (B.21))  for  both  gn  and  gm  in  self-adjoint  form.  We  observe 

(Sm  ~  Sn)  X  9ngm  =  9n  [x  9m)  ~  gm  [X  gn) 

—  \x  ( 9n9m  —  9m9n)\  ■ 
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Tabic  B.l:  Sn  for  Various  Values  of  B 


n\B 

0.2 

1 

5 

10 

1 

0.7593 

1.5708 

2.5704 

2.8363 

2 

4.5379 

4.7124 

5.3540 

5.7172 

3 

7.7511 

7.8540 

8.3029 

8.6587 

4 

10.9225 

10.9956 

11.3348 

11.6532 

5 

14.0804 

14.1372 

14.4080 

14.6869 

6 

17.2324 

17.2788 

17.5034 

17.7481 

7 

20.3811 

20.4204 

20.6120 

20.8282 

8 

23.5280 

23.5619 

23.7289 

23.9218 

C/5 

<D 

05 

> 

CD 


Figure  B.2:  S  for  n  =  1,  2,  3  as  B  Increases  From  1  to  10 


Integrating  leaves 


(sm  -  4)  /  x2gngm dx  =  x2  (gng'm  -  gmg'n)  | 

Jo 

=  9n{l)g'm{l)  -  gm{l)g'n{l) 
=  0 
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which  follows  from  Equations  (B.22)  and  (B.23).  So  for  m  ^  n,  ( gn,9m )  =  0.  To 
normalize  gn  we  choose  Kn  such  that 


1  =  /  x2g2ndx 


=  /  K2sm2(snx)dx 


and  by  trigonometric  identities, 


Integrating  implies 


K 2  f1 

1  =  — -  /  1  —  cos(2sna;)da;. 

^  ./n 


K2 
1  =  thi 

2 

K2 

_  ^ n 

2 


1 - sin(2sr 

2  sn  1 


1 

Sn 


1 - cos(sn)  sin(sn) 


(B.19) 


Equation  (B.17)  allows 


sms,  = 


1  -  B 


COS  Sr 


which  lets  Equation  (B.19)  become 


K2 

n  1 1  _ 


Kl 

2(1  ~B) 


cos 2(sn) 

1  -  B 

[sin2(sn) 


B]. 


Therefore,  let 


AT  = 


2(1  -B) 
sin2(s„)  —  B 


and  for  this  value  of  Kn,  (gn,  gm)  =  bmn-  Notice  that  by  applying  Equation  (B.17), 
linie^!  Kn  =  y/2. 
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We  next  look  for  a  solution  in  the  form  of  the  expansion 


OO 

v(x,  t)  =  7  +  hn(r)gn(x) 

n=  1 

(B.20) 

where  gn(x)  satisfies 

1  r 

~2  x29n(x )  +  S2ngn(x)  =  0 

X  L 

(B.21) 

4(°)  =  ° 

(B.22) 

4(!)  +-B^n(l)  =  0. 

(B.23) 

Further,  with  this  choice  of  v,  the  boundary  conditions  are  satisfied.  Introducing 
Equation  (B.20)  back  into  Equation  (B.7)  gives 

X  X  X  ^  , 

h'n(T)9n(x)  =  -Pi  ~^2(3hn(r)gn(x)  +  ^  Mr)^  x^9n{x ) 

n—  1  n= 1  n—  1 

which  leads  to 

OO 

[h'niT)  +  fan{r)  +  s2nhn(r )]  gn(x)  =  -fa 

n=  1 

and  by  orthogonality  of  {gn} 

K{t )  +  c onhn(r)  =  -  fa  (1,  gn)  (B.24) 

where 

x>n  —  P  +  4-  (B.25) 

Finally,  the  initial  conditions  are  determined  from 

OO 

7  +  K{0)gn(x)  =  ai(x)fi(T*)  +  a2(x)f2(T*) 

n=  1 
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and  the  orthogonality  of  gn  allows  us  to  obtain 


hn(Q)  —  Fn{T*)  —  7  (1,  gn)  (B.26) 

where 

Fn(T '*)  =  fi(T*)  {augn)  +  /2(T*)  (a2,gn) . 

Equation  (B.24)  along  with  the  initial  conditions  produces  the  solution  to  the  first- 
order  differential  equation: 

Mr)  =  [F„(7")  -  7  <!,»„)]  e-“"r  +  f  [-/37  <1,9„)1  d£ 

Jo 

=  Fn(T*)e~u}nT  -  7  (l,gn)  e-WnT  -  [l  -  e“w"Tl  .  (B.27) 

By  utilizing  Equation  (B.25),  Equation  (B.27)  is  reduced  to 

hn(r)  =  FJT*)e-“nT  -  lillM  Ip  +  s2ne~UnT]  .  (B.28) 

UJn 

Next  we  observe  that 

OO 

7  =  'Yjl  0-,9n)  9n{x), 

n— 1 

which  indicates  that  Equation  (B.20)  can  be  converted  to 

OO 

v(x,  r)  =  ^2  [Mr)  +  7  (1, 9n ) ]  9 nix') .  (B.29) 

n=l 

Equation  (B.28)  implies 

K{r)  +  7  <1,  <?n>  =  Fn(T*)e-^T  +  7  (l,gn)  [l  -  ^  -  ^e~^T 

UJn  L0n 

=  Fn{T*)e~UJnT  +  7  (l,gn)  ^  [1  -  e~u}nT~\  (B.30) 


B.10 


where  we  again  make  use  of  Equation  (B.25).  Equations  (B.29)  and  (B.30)  allow  us 
to  state 


v(x,  r) 


OO  OO  2 

J2Fn(r)e—Tgn(x)  (1,  j„>  ^  [1  -  e““q  gn(x). 

Un 


(B.31) 


To  determine  the  average  core  water  concentration  we  find  the  total  core  water  and 
divide  by  the  core  volume  ( Vcore  =  |vrr^) 


r*2n  pit 


Vr 


w(r,  t)r 2  sin(0)d0d6)dr 


core  JO  JO 

3  rrc 


—  /  w(r,  t)r2dr. 


1 0 


(B.32) 


Now  define  xc  =  rc/r0  and  use  the  dimensionless  variables  in  Equation  (B.4)  such 
that  Equation  (B.32)  becomes 


[H20]c(fcr)  =  [H20]m  f  v(x,  r)x2dx. 

xc  Jo 

Through  use  of  Equation  (B.31)  we  obtain 


[H20]c  (tcr)  =  4  [H20]„  {  V  F„(T-)G„(x0)e-“"- 

X"  l  11=1 

00  2 

+7^Gn(l)Gn(zc)—  [1  -  e-— ] 


n= 1 


where 

Gn(x)=  [xegn(m 

Jo 

and  we  observe  Gn(  1)  =  (1  ,gn). 


B.ll 
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